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Abstract 



We study a large class of one dimensional reaction diffusion mod- 
els with quenched disorder using a real space renormalization group 
method (RSRG) which yields exact results at large time. Particles 
(e.g. of several species) undergo diffusion with random local bias (Sinai 
model) and may react upon meeting. We obtain a detailed description 
of the asymptotic states (i.e attractive fixed points of the RSRG), such 
as the large time decay of the density of each specie, their associated 
universal amplitudes, and the spatial distribution of particles. We 
also derive the spectrum of non trivial exponents which characterize 
the convergence towards the asymptotic states. For reactions which 
lead to several possible asymptotic states separated by unstable fixed 
points, we analyze the dynamical phase diagram and obtain the criti- 
cal exponents characterizing the transitions. We also obtain a detailed 
characterization of the persistence properties for single particles as well 
as more complex patterns. We compute the decay exponents for the 
probability of no crossing of a given point by, respectively, the single 
particle trajectories (0) or the thermally averaged packets (0). The 
generalized persistence exponents associated to n crossings are also 
obtained. Specifying to the process A + A — ► or A with probabilities 
(r, 1 — r), we compute exactly the exponents 5(r) and ip(r) characteriz- 
ing the survival up to time t of a domain without any merging or with 
mergings respectively, and the exponents 5a (r) and Vm(?") characteriz- 
ing the survival up to time t of a particle A without any coalescence or 
with coalescences respectively. 9, tp and S obey hypergeometric equa- 
tions and are numerically surprisingly close to pure system exponents 
(though associated to a completely different diffusion length). The ef- 
fect of additional disorder in the reaction rates, as well as some open 
questions, are also discussed. 
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I. INTRODUCTION 



A. Overview 



Reaction diffusion processes are of wide interest in physics, chemistry and biology 
[|TJ. In physics they present a relatively simpler case of non equilibrium stochastic 
processes with non trivial behaviour. Traditionally they have been studied via mean 
field type methods (e.g. law of mass action, local chemical kinetics) M, However, in 
sufficiently low spatial dimension, particle density fluctuations become dominant and 
mean field methods become invalid ||. The role of fluctuations in these processes 
has thus been studied for a while, but has received renewed attention recently as 
new exact results in one dimension Jl[] and systematic renormalisation group studies 
have appeared ||. One interest of these models is their relation to phase ordering 
kinetics via the "coarsening" of domain structures evolving towards equilibrium 
0. In some cases, these can be seen as reaction diffusion processes for defects, for 
instance domain walls in one dimension or XY type vortices in two dimension, which 
diffuse and can annihilate or coalesce upon meeting. These coarsening processes have 
also been much studied recently, especially in an effort to characterize their so-called 
persistence (or survival or first passage) properties for single spins, domains or global 
magnetization H7|-pT[] 

Although many results are now available for reaction diffusion processes in homo- 
geneous situations, comparatively little is known on their dynamics in the presence 
of quenched disorder, which is expected to play a role in many physical realisations. 
It can be introduced in the models in several ways, e.g. in the reaction rates or in the 
single particle diffusion. One can expect that it will strongly modify the behaviour 
of the system in some cases by amplifying the role of spatial density fluctuations. 
These effects are interesting, but difficult to study analytically because of the present 
lack of methods, beyond mean field approximations or perturbation theory, to treat 
the dynamics of such disordered systems. 

Even in the absence of quenched disorder, there is an apparently unlimited va- 
riety of behaviours in reaction diffusion systems. The more complex ones, such 
as oscillatory or chaotic behaviours, become possible for a large enough number of 
species |jl2|-|T5f . In simpler cases, attempts have been made to identify possible uni- 
versality classes and a wide class of models with finite reaction rates, amenable to 
field theoretical treatments, has been studied [|T^,|T^]. For instance, branching and 
annihilating random walks (BARW), i.e reactions such as A — > mA and A + A — > 
or A — > exhibit transitions from inactive (no particle) to active states, which 
were found to be either in the universality class of directed percolation [p~8| pOf| (odd 



number of offsprings) or in the so called parity conserving class (even number of 
offsprings) ]T7|] . This was confirmed by exact results in one dimension |H[ . Related 
types of models describe epidemy propagation, such as A + B — > IB (with rate k) 
and either B — > A (recovery) or B — > C (immune) (rate 1/r) were also studied via 



RG [E^] (see [23[ for review). The effect of quenched disorder has been studied in 



this class of BARW models, via random rates k(x) and t(x) but with limited success 



as the RG flows to strong coupling p3| . As for directed percolation with disorder, 
it is still a largely open problem [^-^ . 

There is a simpler class of homogeneous models without branching (i.e without 
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particle production), such asA + A^O, A + A^A J27|,[28|], A + B -> [||.[3C 
etc.. which has still non trivial behaviour [31]. One interesting phenomenon is 



that in low enough dimension, the process becomes diffusion limited rather than 
reaction limited. Indeed particles in close proximity react quickly and the remaining 
particles are typically separated by a length related to the pure diffusion length 
Lo(t) ~ (Dt) 1 / 2 . This leads to a decay of specie density, e.g. ri^(t) ~ t~ d l 2 for 
d < 2 in the case of A + A — > 0, slower than the mean field prediction ~ t^ 1 

valid for d > 2 (for A + B — > a related argument leads n^(t) ~ t~ d ^ for d < 4) 
This type of results for such models in the pure case are well established from 



heuristic arguments, numerical simulations, perturbative RG |27|j28| , |32[| and in some 
cases rigorous methods |33| , |3ll . It is now interesting to investigate how disorder will 
modify these behaviours. With disorder, models in this class are easier to study 
than the BARW type models, although it is still a difficult task. The reactions 
A + A — > and A + B — > have been studied using perturbative field theoretic RG 
methods for particles diffusing in random flows, either in two dimension p4] , |35| or 
in a special hydrodymamic flow ||36||. As can be expected from the study of single 



particle diffusion in such flows p7| , |38| , the behaviour should be qualitatively different 
in the case of potential disorder, which tends to segregate the particles and slow 
the reaction (and the diffusion) than for hydrodynamic flows which tend to mix the 
particles and increase the effective reaction rate (and leads to hyperdiff usion) . The 
competition which arises when both components are present has been studied very 
recently in d = 2 in [34] and in d = 2 — e | 39| . Remarkably, the one dimensional 



problem seems quite far from the reach of such perturbative RG methods and no 
generic result is available at present |^(| in that case, hence the interest of the 
present study. Indeed, in d — 1 only potential disorder can exist and is known to 
lead to ultraslow single particle diffusion described by a strong disorder (i.e zero 
temperature) fixed point [41] . To make progress in d = 1 requires developing non 
perturbative techniques, which is the aim of the present work. 



B. Model 



In this paper we study a broad class of reaction diffusion models where particles 
diffuse on a one dimensional lattice and can react or annihilate upon meeting. Apart 
from their reactions the particles are non interacting. More specifically, each site 
of the lattice can be in one of several possible "states", labeled k = 0,1,. .n — 1. 
k = corresponds to the empty state with no particle present at that site, k = 
l,..,n — 1 correspond to the presence of particles of different types. When two 
particles (i.e states) k\ > and k 2 > meet they react and give a another state 
k with a probability k . k may be the empty state k = 0, corresponding to an 
annihilation. The reaction is thus a stochastic process 



k\ + k 2 — > k with probability 



Km 



characterized by a fixed transition probability matrix which satisfies 



(2) 
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The matrix W kl k2 can be extended to include k\ = by defining 

Wj* = W k k , fl = 5 k>k , (3) 

for any k, k', which is the property expected for an empty state (A + — > A with 
probability 1). 

One prominent example will be identical particles A which react upon meeting 

as 

A + A — > with probability r 

A + A — > A with probability 1 - r (4) 

In this case there are only two states: k = corresponds to no particule present (0) 
and k = 1 to one particule present A. The transition matrix is then W® = W 7 ^ x = 

W?o = Wjo = 0, Wji = Wjo = 1 and = r, = 1 - r. 

We will obtain results for processes within the above class (jl]) and study some 
specific examples. We will restrict to symmetric reaction rates Wt k2 = Wt k . 
Asymmetric rates, depending on the side from which the two species come in contact, 
can be defined in d = 1 and can be studied by the same methods. We will mostly 
consider reaction diffusion processes with only a finite number of possible states (or 
species). Processes with unbounded number of states (n = oo) can also be studied 
by the present method, and we will give some examples. Classifying the much larger 
variety of complex behavior possible in that case is beyond the scope of the present 
study. Other extensions include randomness in the reaction rates, which we will 
briefly discuss in the end. 

Up to now we have not specified the way in which particle diffuse, nor the reaction 
rates. Let us first concentrate on the process (f|) and recall the known results in the 
case of pure diffusion (i.e homogeneous hopping rates), which has been extensively 
studied. It is of particular interest in one dimension since it is also a model for 
zero temperature domain growth in the ferromagnetic g-states Potts model (with 
Glauber dynamics) where r = l/(q — 1) ||11|| . The case q = 2 (Ising) corresponds 
to walkers (i.e domain walls) always annihilating when they meet and q = oo to 
walkers always coagulating ||42|| . It is known that the reaction rate can be chosen 
infinite (immediate reaction upon meeting) without changing the universality class, 
and the same will hold here in presence of disorder, hence our general choice of 
model ([]]). For all r the concentration of particles A is known to decrease as: 

n A (t) ^ c(r)(Dt)- 1/2 (5) 

where the (r — dependent) coefficient is expected to be universal (e.g c(l) = (87r) -1 / 2 , 
[pp. More detailed properties, such as persistence, have also been studied. The 
probability S(t) that no particle A (domain wall ) has crossed a given point O up 
to time t has been shown to decay as: 

S(t) ~ L (t)- d ^ (6) 

where Lo(t) ~ \fi is the characteristic length and 9(r) the so-called persistence 
exponent. S(t) also corresponds to the probability that a spin has never been flipped 
up to time t in the Potts model. The exact expression of 9(r) was obtained in [0,0. 
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The domain size distribution has also been computed for this process in [45]. Two 
new independent exponents, ip and 5 were introduced and studied in [0,0 to 
characterize the persistence (survival) of domains for this model, but up to now no 
exact result is available for these exponents. The concept of persistence properties 
was extended to other observables, finite temperature and studied in a variety of 
other models: persistence for global order parameter in [jHj, spin block persistence 
in [58], generalized persistence and large deviations in P9"|JoT)|] and persistence for 



fluctuating interfaces [51 



In this paper we study the case where the hopping rates are inhomogeneous with 
short range correlations, corresponding to random local bias. The generic model 
for this type of disorder is the Sinai model where each particle performs Arrhenius 
diffusion in the same energy landscape U n where the local random forces U n — U n+ \ 
are independent random variables, of zero average (we restrict the present analysis 
to zero global bias). Various analytical results are known for the single particle 
Sinai model |52|-p8|. Diffusion is ultraslow as x ~ L(t) = (hit) 2 . Recently we have 
reexamined this model 



59,60 



using a real space renormalization group method 
(RSRG) which yields exact results at long time. In the present paper we apply the 
RSRG method to study reaction diffusion of the type ([!]) for particles in a Sinai 
landscape. Some of the results have already appeared in [59]. Although we give 
here a detailed treatment of the reaction diffusion RSRG, we will rely on Ref. |Hj 
for all details concerning the single particle diffusion aspects of the problem (which 
we will only sketch, referring the reader to |Kj for details). Note that we consider 
here only models where all particles share the same diffusion property (i.e. see the 
same landscape, and have same diffusion coefficient). Thus this does not include 
reaction diffusion models such as Ising domain walls in a random field, for which a 
specific treatment is necessary and which are studied in ]^9| , |6T| . Similarly, relations 
to other problems such as disordered quantum spin chains p^j-po] or disordered free 
fermion models are discussed in ]^,Q. In particular we have chosen to discuss our 
present results exclusively in terms of reaction diffusion dynamics, and not in their 
equivalent formulation as (non hermitian) ID quantum models (see |66| for details 



of such relations in the pure case). Finally, note that an exact RG has also been 
applied to the problem of coarsening of the pure ID $ 4 model at zero temperature 
for which persistence exponents have been computed f7],|7j. 

As for the single particle problem, the RSRG method allows us to compute a 
number of quantities, and, remarkably, even some which are not known for the 
corresponding pure model (e.g. the domain persistence exponents 8 and ip). This 
makes the disordered case all the more interesting to study. We find that reaction 
diffusion processes in a Sinai landscape are strongly controlled by the ultraslow 
diffusion, e.g the relevant length scale is the diffusion length L(t) ~ (hit) 2 , but 
that they still possess non trivial reaction properties. We will characterize a broad 
set of universality classes, containing all reactions of type (|I]). The reaction times 
(provided they are finite) do not affect any universal quantity, so that we can consider 
the reactions as instantaneous for practical purpose. As discussed in |Kj there 
are other single particle diffusion models with short range correlated disorder in 
one dimension apart from Sinai's model universality class, such as random barriers 
(symmetric hopping rates) or random wells. For interesting behaviour to occur, 
however, algebraically broad distributions are required from the start. Some results 
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for reaction diffusion processes with this type of single particle diffusion have been 
obtained in [|68|| . 

The outline of the paper is as follows. In Section |I| we detail the RSRG method, 
first recalling known results in the case of a single particle in [II A| , then deriving the 
RSRG equation for reaction diffusion models [11 B[ The fixed points of this equation, 
and some physical properties of the corresponding asymptotic states are studied 
in |11 C] and |11 JJ| respectively Section [TIT] is devoted to a detailed analysis of the 
dynamics near attractive or repulsive fixed points and of the convergence towards 
the asymptotic states. Throughout the paper we apply our results to the process 
(fD but in Section [TV] we discuss some applications to other examples of processes. 
In Section |V| we study persistence properties. Section [VI I| contains the conclusion. 
Some more technical but useful details are contained in the appendices. 



II. RSRG METHOD FOR REACTION DIFFUSION AND ASYMPTOTIC 

STATES 

A. RSRG for Sinai landscape and single particle diffusion 



The model for the diffusion of a single particle in one dimension can be defined, 
with no loss of generality |5^| , as the Arrhenius diffusion in a "zigzag" potential U(x) 
represented in Fig. |I] (a). It consists in a set of bonds, each bond (between Xi and 
Xi + \) being characterized by an energy barrier Fi = \Ui — i7i+i| (where U{ = U(xi)) 
and a length U = \xi+i — Xi\. The energy landscape is chosen by choosing a pair of 
bond variables F, I, independently from bond to bond, from a distribution P(F, I) 
normalized to unity. 




1 1 2 1 3 

FIG. 1. (a) Energy landscape in the Sinai model (b) decimation method: the bond 
with the smallest b arrier F m { n — F 2 is eliminated, resulting in three bonds being grouped 
into a single one. 



The RG procedure, which captures the long time behaviour in a given energy 
landscape, is illustrated in Fig. [I] (b) and consists in the iterative decimation of 
the bond with the smallest barrier ||59|| , say F2, and to replace the three bonds 
1,2,3 by a single renormalized bond with barrier F' = Fi — F 2 + F 3 and length 
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I' = l\ + I2 + I3. The new variables remain independent from bond to bond. To 
write the corresponding RG equation it is convenient to introduce T as the smallest 
remaining barrier at a given stage of the decimation and the rescaled variables 
r] = (F — r)/r and A = l/T 2 . The RG equation for the probability distribution p9| 
P r (?7, A) reads 

(Td r - (1 + r))d v - 2Xd x - 3)P r (r/, A) = P r (0, .) * A P r (., .) *„ >A P r (., .) (7) 



and coincides with the one derived in [^3J for the closely related problem of disordered 
quantum spin chains. The symbol * A denotes a convolution with respect to A only 
and with respect to both rj and A. The probability distribution is normalized 
to unity as drj J +o ° d\P r (r], A) = 1. 

The landscape is characterized by the large scale variance of the potential: 



([/. _ u.y ~ 2a\k^\ (8) 

with the distance from site i to site j, which is exactly preserved by the RG. 
Thus we will set a = 1 in the following. Restoring a simply amounts to a rescaling 
of lengths, and in particular a drops out of all (universal) ratios of lengths that we 
consider later. As shown in |63|j6^1 , the RG equation (0) leads at large T (using 
Laplace transforms) to the following fixed point P*(t], A) : 

P*( V , A) = LT;\ (-^= e-^-h ^ (Q) 

Thus for large V one finds that the average bond length Ir and the number of bonds 
n? per unit length are |7D| , respectively: 

h = ^ 2 n r = ^ (10) 

The renormalized landscape allows to study the dynamics of a single walker 
starting from a given point O at t — 0. The decimation of barriers smaller than 

r = Tlnt (11) 

corresponds to the elimination of (logarithmic) time scales shorter than the Arrhe- 
nius time t for the particle to cross the barrier. We are choosing everywhere time 
units such that the (non universal) microscopic attempt time scale to be set to unity 



(arbitrary units can be recovered by setting T = T\n(t/to) in what follows ||71|| ). 
Since at long time (i.e large T) the renormalized landscape consists entirely of deep 
valleys separated by high barriers, a good approximation to the long time dynamics 
is to place the walker at the bottom of the renormalized valley at scale T = Tint 
which contains the starting point O, since with high probability it will be near to 



that point |52]. Upon proper rescaling of space and time this approximation be- 



comes in fact exact as T tends to +00. This defines what we will call the "effective 
dynamics" in the following and is illustrated in Fig. [I]. This allows to recover the 
scaling x ~ (lnt) 2 for the single particle diffusion as well as many other exact re- 



sults detailed in |59|,|60[|. Since it is customary, when studying reaction diffusion 



processes, to compare densities of reactants with a characteristic scale of diffusion, 
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we give here the the exact expression for the single particle r.m.s displacement, or 
"diffusion length" at large time: 



<* 2 (t)>^W^T 2 (lnt) 2 (12) 
o V 5 

To study reaction diffusion processes it will be necessary to consider "valleys" 
(two consecutive bonds sharing a common potential minimum). We thus slightly 
generalize the above RG equation (|?D to follow the distribution of renormalized 
valleys. The RG equation for the valley probability distribution -P r (771,772) a t scale 
T in rescaled variables (7/1,772) reads 

TdrP r ( Vl , V2 ) = [(l + Vi)d Vl + (1 + 7/2HJ +2] P T (7/1, 7/2) 

+P r ( Vl , .) *, 2 P r (0, .) + P r (., 0) * m P r (., m ) (13) 

where we have omitted (i.e integrated over) the lengths for simplicity. The large 
time (r) behaviour of this equation can be studied similarly. Valley distributions 
which have the decoupled form P r (7/1, 7/2) = P r (i]i)P r (rj2) where P r (7/) satisfies the 
bond RG equation (0) are of course solution of the RG equation for valleys fll3"l) . The 
subspace of such decoupled distributions (called decoupled subspace in the following) 
is thus preserved by RG. Since the initial condition is uncorrelated, the RG flow 
defined by (0) remains in this decoupled subspace, and converges towards the fixed 
point P*(r]i,i]2) = P*{f]i)P*{ri2) = e~ m ~ m . This convergence result extends to the 
case of small correlations between valley sides as will be discussed below. 

It was shown in |M that the convergence towards the bond fixed point P* (77) = 
e~ v is like ^ with eigenvector (1 — r])e~' n . Thus the convergence towards the the 
valley fixed point P*(rji, r] 2 ) = e~ vi ~' n2 within the decoupled subspace is of the form: 

(e^ 1 + £(1 - m)e- m + -.) (e"" 2 + £(1 - m)e~ m + •••) = e~"(l + £(2 - Vl - t/ 2 )) + ... (14; 

and is also as f with eigenvector (2 — 7/1 — r\ 2 )e~ r]x ~ m . 



B. RG equations for reaction-diffusion 

We now turn to diffusion reaction models of type ([!]) in one dimensional land- 
scapes with random local biases. 

From the results on the dynamics of a single particle in Sinai landscape recalled 
in the preceding section, it is clear that one can study most of the properties of 
the initial reaction diffusion problem by following its evolution under the effective 
dynamics. It also becomes obvious that one must now consider valleys, and the 
species contained in these valleys. At the decimation time scale T = Tint in some 
places in the system, two valleys will merge into one and the reaction (|TJ) governed 
by the rates matrix W will take place. This process is illustrated in Fig. |2]. The 
errors made by this approximation are expected to become again smaller at large 
time, as will be discussed later on. 

The general method to study the process ([[]) is thus to associate to each valley 
the specie which it contains, which is one of several possible states k (k = being the 



FIG. 2. (a) macroscopic state of the reaction diffusion process: each renormalized 
valley is either empty (k = 0) or contains a particle of type k > (b) at time scale t 
such that the barrier F = T = Tint is decimated, the state (particle) k\ in the decimated 
valley moves to the neighboring valley and reacts with k 2 to produce k with probability 
Wt , , as the two valleys are merged into a single renormalized one containing k 



empty state). A convenient initial model is thus one where each valley and its content 
is statistically independent and characterized by a probability distribution P^ivii V2) 
with J2k Pk(ViiV2) — P T (Vi: Vi)- It remains so under the RG. The effective dynamics 
is described by the RG equation: 

Td r P k r ( Vl , r, 2 ) = ((1 + m )d Vl + (1 + m)d m ) + 2)P fc r fa 1 , r] 2 ) 
+Wt lM [P£ l (m, •) ** PkM •) + °) *m P U;V2)] (15) 

where summation over repeated indices is implied. The summation over k yields 
back the valley RG equation ([13]). Since the average length of a valley is 2lp, The 
total concentration n&(t) of a given specie k (the total number of fc-particles per 
unit of length) is given as: 

Mt) = \ n rPl (16) 

where: 

pl= [ Pk(m,V2) (17) 

and nr is the number of remaining bonds (|ToD per unit length at scale T = Tint. 
We stress that the RG equation flTSP is more complicated to analyze than ( JTB| ) since 
it cannot in general be factorized into bond distributions. 

However, it turns out that there is still a simple subspace of distributions which 
is exactly preserved by the RG equation ( Jl5| ) . It is the subspace of functions of the 
sum r/i + r\2. 

P^r ]l ,r ]2 )=Hl(r ] = r ]l + V2 ) (18) 
where the functions Hk(rj) satisfies: 
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Td r H T k ( V ) = ((2 + V )d v + 2)Hl( V ) + Wl M Hl{.) H T k2 {.) (19) 

which conserved the normalization ^ k p\ = f^dr] vJ2k^k(v)- This subspace 
plays an important role in the following. Already one sees that both the fixed point 
e -m-V2 anc [ kh e leading eigenvector (2 — 771 — ?]2)e _r ' 1 ~ 7?2 of the linearized landscape 
RG equation fll3|) belongs to this subspace. 



C. Fixed points of the RG equations and asymptotic states 

We now determine the fixed point solutions of the RG equation for the valley 
distributions (|15|). We already know that the sum ^2 k P k (rji, r/ 2 ) = -P r (771,772) con- 
verges towards the fixed point of equation ( p~3| ) P*(r)i, i] 2 ) = e~^ r,1+m > which describes 



the landscape. It is thus natural to look for fixed points of (|15|) of the following form 

Pkfaih) =Pl P*(vnV2) =ple-<*+ri (20) 

where p* k > and ^2 k pl = 1 by normalization. Plugging this form into flip]) leads 
to a consistent T independent solution if the pi satisfy the condition: 

Pi = W k k lM pl lP l 2 (21) 

Note that any solution of this equation satisfies a priori J2kP*k = or 1, as a 
consequence of (|3|). Thus apart from the unphysical solution of fl2T|) where all p* k 
vanish, all other solutions are automatically correctly normalized. 

In general the equation (^TJ) has several solutions and thus there are several fixed 
points to the valley equations (|T5|). Clearly some of these fixed points are attractive 
and correspond to possible large time asymptotic states for the reaction diffusion 
process while other fixed points are repulsive. In some cases several attractive fixed 
points can coexist and lead to a non trivial phase diagram. 

The stability of each fixed point, as well as their convergence property, will be 
studied in detail in the next Section. Here we just mention one important result. 
The dynamics in the vicinity of a fixed point pi is determined by the following 
stability matrix: 

M kj = Wljfr (22) 

We denote by /i a the eigenvalues of M, and by the associated eigenvectors 
Mkjp'j = Ha.pt- F° r n—1 reacting species (in addition to the empty state k = 0), M 
is an n-dimensional matrix with non-negative elements. One eigenvalue is trivial, 
being simply \i\ = 1 of eigenvector p k =1 oc p* k (|2l|). The other n — 1 eigenvalues // Q 
(a = 2, ...n), which can be complex in general, have a smaller modulus, due to the 
Perron-Froebenius theorem. We will focus in the following on simple enough pro- 
cesses where all real, but we will also give example of more general reaction 
diffusion where complex eigenvalues arise (e.g. reactions with cycles). 

Let us consider here only fixed points with all fi a being real. The result of the 
next Section is that, in that case, a given fixed point if attractive is the n — 1 
eigenvalues fj, a < 1/2 for a — 2, ...n. It is repulsive if at least one of these /i Q > 1/2 
(and is then repulsive along the corresponding eigendirection). 
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For example, in the case of the reaction-diffusion process (f|), there are two 
solutions of ( pTj) and thus two fixed points. One is the empty state E=(j»q = 1 
, p* A — 0) and the other is S=(pg = r/(l + r), = 1/(1 + r)). The stability 
matrix associated to E is simply the 2x2 identity matrix (i.e. fii = 1 = /i2) 
and this fixed point is thus repulsive. The matrix associated to the state S reads 
M s = ((^i, ^x), (^x, witn eigenvalues /xi = 1 and /j 2 = 0. The fixed point 
S is thus attractive and corresponds to the asymptotic state which represents the 
large time behaviour of the system. 

In fact the reaction-diffusion process (|j) possesses an interesting property: the 
outcome of a sequence of reactions does not depend on the order it was performed. 
We call these processes "associative processes" . They have the property that M 2 = 
M, i.e /i = 0, 1. Some properties of these associative processes are detailed in 
Appendix (||). 



D. physical properties of the asymptotic states 



We now study some physical consequences. Each attractive fixed point corre- 
sponds to a possible large time behaviour of the system, i.e an asymptotic state. If 
there are several attractive fixed point, the one chosen by the system will depend 
on the initial value of the parameters (mainly the specie concentrations). 

From the results (^) and (|T0| ) we obtain that in an asymptotic state (charac- 
terized by a set of p* k solution of (|21~D , the density of specie k behaves at large time 



as: 



Pi 



THn\t/t c 



(23) 



where we have restored the microscopic attempt time scale to |72|- Note that this 
result (|23"D represents the leading large time contribution, subleading corrections 
(which become dominant only if pi = 0) will be determined in the next section. 
Interestingly, this leading behaviour is independent of the initial concentration (pro- 
vided it is in the basin of attraction of the fixed point). This universality property 
can be further characterized by computing universal amplitudes. In pure models a 
commonly studied amplitude is the product of specie concentration by the diffusion 
volume. In a disordered model one has more choices of definitions, but we will define 
an amplitude as in [|39|. Here we find the exact result for the following universal 
amplitude (associated to specie k): 



(24) 



This gives for instance A(r) = 1/(1 + r) | y y for the process 
for A + A —>■ studied in [BO] by perturbative methods [[T3| . 



i.e A = 0.291071.. 



From the statistical independence of valleys in an asymptotic state, information 
about the spatial distribution of the species can also be obtained. For instance one 
can define "domains" , in the simplest intervals between particles (i.e between 

non empty states k ^ 0, irrespective of their content) We can now compute exactly 
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the distribution of size of "domains". From the above form of the fixed points the 
normalized distribution D-p(l) of domain sizes / in an asymptotic state takes a scaling 
form Dr(l) = ^D*(X = l/T 2 ) which can be computed as follows. 

The above RG equation for valley distributions (|15D can be readily extended to 
Pk(vi-> ??2, Ai, A2) which takes into account the rescaled lengths Ai = h/T 2 , X 2 = h/F 2 
of the two bonds of the valley, extending (|7|). The generalized fixed point (p0|) reads: 

P:( Vl , V2 ,\ 1 ,\ 2 ) = p* k P*( Vl ,\i)P*(V2,\2) (25) 

where P*(rj, A) is the fixed point solution (§) of the bond RG equation. A domain as 
defined above, is thus a set of consecutive empty valleys between two occupied val- 
leys, together with one bond in each of the occupied valleys (see Fig. ||). Since in the 
asymptotic state valleys are statistically independent and either empty (k = 0) with 
probability Pq-P*(Ai)-P*(A 2 ) (where P*(X) = J P*(X)) or contains a particle (k ^ 0) 
with probability (1 — £>q)P*(Ai)P*(A 2 ), one easily obtains the Laplace transform of 
D*(X) as: 

= f^rrm = ji^g = ^f^ i (26) 

where we have used the explicit form @ for the fixed point bond distribution. 



Formula (26) can be inverted and yields the distribution of rescaled domain sizes: 

+00 

D*.(A) = tana ^ (a + mr) e - x{a+nn)2 (27) 



71= — OO 



tana m 2 
= 2 , msin(2«m)e * (28) 

* m=— 00 

with a = arccos^Q. 

This can be applied, e.g. in the case of the process (f|). Substituting^ = r/(l+r) 
in (p8|) yields the distribution of distances between neighboring walkers A. Note that 
the case where walkers A always coalesce upon meeting (r = 0, po = 0) corresponds 
to a — > I and in this limit ( p8|) becomes: 

+00 

D;, q=0 (X) = (2Avr 2 (n + 1/2) 2 - i) e -^ 2 (i/2+n)2 (2g) 

n=— 00 

o +°° 2 

^ (_l)"+i m V^ (30) 



0FA 3 / 



2 

m=— 00 



It is interesting to compare the result ( p8|) concerning the disordered case with 
the result of Derrida and Zeitak p5[ for the case of homogeneous hopping rates. For 



small domain sizes (A — > 0), the distribution vanishes much faster in the disordered 
case (as ~ A -3 / 2 exp(— 1/A)) than in the pure case ( as ~ A). For large domain 
sizes (A — > 00), both have exponentially decaying behaviour (except for q = +00, 
i.e r = in the pure case). In addition, in the present case the consecutive domains 
lengths are statistically independent, which is not the case for the pure system. 

The above calculation is easily generalized to compute the distribution of relative 
distances between two walkers of a given species k, simply by substituting p* — > I— pi 
in the above formula fESI). 
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III. DYNAMICS NEAR FIXED POINTS AND ASYMPTOTIC STATES 



In this Section we study the dynamics near the possible fixed points of the valley 
RG equation (attractive and repulsive). 

We will first focus strictly on the effective dynamics exactly described by the RG 
equation (flq), and mention some possible corrections in the real dynamics at the 
end of the section. 

For the effective dynamics we will solve the problem in two steps. As mentionned 
above the matrix M in ( |22| ) and its eigenvalues fi a control the asymptotic dynamics. 
Interestingly they readily provide an approximation of the dynamics, which we will 
call the "uniform approximation", which is interesting as it allows to classify the 
spectrum of eigenperturbations and, in the case of real eigenvalues, already allows 
to see whether a given fixed point is stable or unstable. 

In a second step we will obtain the exact results for the spectrum of eigenper- 
turbations. 



A. First step: uniform approximation 

It is natural to define the total occupation probability of specie k at scale T = 
Tint as p\ = f r]ir)2 Pk(Vi,V2)- The difficulty of the problem comes from the fact 
that it does not satisfy a closed equation. However, if one also introduces p^(0) = 
I m -ffcXO) V2) = / 1 PkiVii 0)) i- e t ne occupation probability of specie k of the valleys 
just being decimated at T, one can obtain a closed coupled equation by integration 
of (|T5| ) which reads: 

Tdrpl = 2(-p£(0) + (31) 

It is then tempting to set, as an approximation, p£(0) = p\. This would be correct 
at any of the fixed points (p0| ), but since we are studying convergence to a fixed 
point, it is an approximation which amounts to neglect the dynamical correlations 
between the deviations in the specie concentration and in the distribution of barriers 
heigths. For this reason we call it the "uniform approximation". It yields the 
following approximate closed RG equation for the p\: 

Tdrpl = 2(-p£ + W£ xM pl x pl 2 ) (32) 

which preserves the normalization condition Y2Pk — 1- This approximate flow has 
the same fixed points p\ = p* k than the true one This equation, remarkably, 

is reminiscent of a "mean field type" rate equation, except that the role of "time" 
would be played by the variable ln(Tlnt). 

The relaxation of (|32|) towards any of these fixed points is studied by setting 
Pk = P*k + fk an d linearizing for the small perturbation f£ around the fixed point 
p* k . It yields, in terms of the matrix M introduced in (|22|) : 

Td T f k =2{-f k +2M kk ,fl,) (33) 

and thus the convergence towards the fixed point has components behaving as r~ Aa 
where the exponents are given in terms of the eigenvalues fi a of the matrix M as 
A a = 2(1 — 2 Ha) with a = 2, ..n. 
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So this "uniform approximation" would indicate that a fixed point is stable if all 
Re(fi a ) < 1/2 for all a = 2, ..n, and unstable otherwise. Remarkably, this stability 
criterion coincides with the exact result when the eigenvalues are real as we will 
now show, even if the naive convergence eigenvalues A a are not correct (they are 
"renormalized" to larger absolute values). 



B. Second step: full dynamics near a fixed point 

Up to now we have studied the convergence of the landscape alone (|13"D, and the 
convergence of the p r within a uniform approximation assuming p r (0) = p r . 

We now study the full dynamics near a fixed point solution of the full reaction- 



diffusion equations (|i~5p. We will indeed find that there are some correlations be- 
tween deviations in total occupation probabilities (from the fixed point concentra- 
tions) and deviations in the barrier distribution profile (from the simple fixed point 
shape e~ v ), resulting in deviations with respect to the uniform approximation. 
We thus consider a perturbation around the fixed point of the form 

P£(Vi,V2) = if k + c r k ( VuV2 ))e-^) (34) 
and linearize the equation for the perturbation (^(771,772) 

r<9 r c[(77i, r/ 2 ) = ((1 + 771)^ - 771 + (1 + 772HJ - m)c r k (vi, V2) (35) 

pi rV'2 

+M klM [ drf((f kl (rf,f j2 ) + (f kl (rf,0))+ / drf{^^rf) + ^rf)\ (36) 



Jo 



where we have used the symmetry of the W, and the definition of the matrix M 

(H). 

Note that in the end we are interested into the behavior of the species proportions 
Vl = I PkiVi, ifr) =Pl + e r k where e r k = [ c T k ( Vu V2 )e~^ + ^ (37) 

J m ,m J ni ,m 



The normalisation condition of course implies that Ylk e \ = 0- 

Decomposing c^(r/!, 772) upon the eigenvectors p£ corresponding to the eigenvalues 
/i Q of the matrix M as c^r/^r^) = c^(r/i, r/ 2 )p k , we obtain decoupled equations 
for the coefficients 0^(771,772) 

rdr&rn, 772) = ((1 + Vi)d Vl - r?! + (1 + V2 )d V2 ) - mKiVi, t/ 2 ) (38) 
H«[ dr j , (c r a ( V ',r l2 ) + c r a (r l , ,0)) + / d V '(c r a ( Vl , 770 + ^(0,77')] (39) 



JO JO 

For a given fi a , we look for solutions behaving as ~ r~* Q and determine the 
exponent $ a as a function of the eigenvalue 71^. Here, a priori both fi a and § a can 
be complex. 

Before we study this equation for general fi we will first study the simpler cases 
fi = and fi = 1. Note that for associative processes this will be sufficient. 
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1. Study for \x a = 

This case is important for naively stable fixed points of associative processes 
(which have all \i a = for a = 2, ...N — 1). The T dependent equation fl3§D with 
H a = can be integrated out explicitly starting from its initial value at T 

cl(vx,V2) = (^) V^-^+^c^la + m) + m) - 1) (40) 

1 1 1 

and thus the convergence of proportions p\ towards p* k is governed by fl3~T|) 

poo poo poo poo 

e v a = / d m / efy 2C r( M2 ) e -fe+^ = / d Vl / dy 2 c r a °(y 1 ,y 2 )e-^ + ^ (41) 



.7 — 1 — 1 



So for generic initial perturbations c^°(?7i, 7/ 2 ) that are not exponentially growing as 
Vi,2 — > 00 ) we obtain that the convergence towards the fixed point is exponential in 
r. Note that here this exact result is very different from the naive approximation 
which would predict a convergence as T" A with A = 2. To understand why this is 
so, one can compute from the exact solution (|41"D the ratio: 



pW)-p1 r 
pI - pI r, 



o 



(42) 



which is found to grow with T! This is why the uniform approximation is particularly 
bad for this case where it predicts a power-law instead of the exponential convergence 
in T. 



2. study for fi a = 1 

To study the dynamics ( p8|) in the case \i a = 1, it is useful to introduce the 
function h^Jj}\,Tfy) = d^d^c^rji, 772), since it satisfies the closed simpler equation 

TdrhKm, 772) = ((1 + Vl )d Vl + (1 + 772KJ + 2 - 77! - 772)^(77!, 772) (43) 

that gives after integration from an initial condition at T 

^(t?!, 772) = (^) 4 e-^- 1) ^ + " 2+2) ^(^(l +771) - 1,^(1 + %) - 1) (44) 
1 1 1 

and thus for initial conditions /i^°(r/i, r/2) that are not exponentially growing at 
771,2 —* 00, we obtain that the the function converges towards exponentially in 
T. This means that the perturbations converges exponentially in T towards the 
decoupled subspace 

c T a (Vi,V2)=ipl(Vi)+^(V2) (45) 

We now study the convergence towards the fixed point in that decoupled subspace 
to see if there are solutions behaving as 7/^(77) ~ Y~^ a ip ol {vi) : the equation for ^(77) 
reads 
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- ®aMv) = [(1 + V)d v - rjlMv) + 2 ^ d V 'Mv') + #a(0) (46) 

Jo 

The only well-behaved solutions are found to be: 

$ Q = +1 ^) = i(0)(l-f|) (47) 

$ a = -2 1>(7 } )=il> a (p)(l + 27 t ) (48) 

The first solution corresponds to the convergence as 1/T of the landscape dis- 
cussed previously. This is not surprising since the linearized RG equation ( |13|) is 



exactly fl3"8"|) with /i = 1. Note en passant that (£3) shows that small correlations 
between barriers in the same valley decrease exponentially fast in T towards the 
subspace of statistically independent distributions. Since this landscape eigenvector 
satisfies e T a = f c^ivii i] 2 )e~^ Vl+ri2 ^ = it does not affect the species proportions 
p\. It is the second eigenvector in fl4*8] ) which is relevant for the reaction diffusion 
processes since ^ 0. It corresponds to the unstable eigenvalue (growth as Y 2 ) 
associated to a naively unstable fixed point (e.g. of an associative process). Note 
that in that case, the unstable eigenvalue found $ = —2 coincides with the naive 
value A = 2(1 — 2/x) = —2 of the uniform approximation. 

Physically, this eigenvalue can be understood for e.g. the process (Q). The 
unstable fixed point is the empty state E with p = 1, p A = 0. Now if one starts at 
t' very close to the fixed point, there are very few A and their number should not 
vary at first, as they will rarely meet. This is indeed exactly what the above results 
says, namely: 

n A {t) = \ v p T A ~ j^T- t iP\ + (paW ~ P*a) 1 ^) (49) 
and using that E has p* A = 0. 



3. Study for general \x 

We now study the case of a general It turns out that one can find solution 
of the original equation (j3£f) under the form 

c£(»7i,7fe) =r-*°H Q (7 il +rj 2 ) (50) 

where the function H a (z) satisfies the differential equation 

= (2 + z)Hl{z) + ($ a + 1 - z)H' a {z) + (2^ Q - l)H a (z) (51) 

together with the boundary condition 2H' a (0) + § a H a (0) = 0. 

The only well-behaved solution at z — > oo is the confluent hypergeometric func- 
tion H a (z) = U(l — 2/i a ,3 + $ Q ,2 + z), and the boundary condition at z = 
determines the possible exponents <3>q, that should satisfy : 2U'(1 — 2fi a , 3 + $ a , 2) + 
$^(l-2/i a ,3 + $ a ,2) = 0. Using the identity zU'(A, B,z) -(z + l-B)U(A,B,z) = 
—U(A — 1, B — 1, z), this equation for $ a reduces to 

U{-2(t a ,2 + $ a ,2) = (52) 

We now discuss the behaviour of the solutions of this equation. One must dis- 
tinguish two cases: 
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4- real \i 



Let us start with \i real. As mentionned above, one must have \i < 1. For 
fj, = 1 the equation (52) reduces to <£> 2 + $ — 2 = which admits the two roots 
$ = 1 and $ = —2 and one thus recovers the eigenvalues (48). The equation (52 ) 
continues to admit two finite roots when \x belongs to the interval | < \i < 1, which 
we denote $ + (/x) and $~(/i) (with <3> + (l) = 1 and $"(1) = -2). The behaviour 
of these roots as a function of \i is plotted in Fig. (||) and ([|). As /i is decreased 
from 1, $ + (l) increases and diverges when \x — - > l/2 + while $~(/i) increases from 
—2 to $~(l/2) = 0. For /i < 1/2 the equation ( |52"D admits only one finite root 
$(/i) = $~(/x) which is positive and with $(1/2) = and $(/i) — > +oo as ji — > + . 



4.0 r 

3.5 - 

3.0 - 

2.5 - 

2.0 - 

1.5 - 

* 1.0 - 

0.5 - 

0.0 - 

-0.5 - 

-1.0 - 

-1.5 - 

-2.0 - 

o.c 

FIG. 3. Plot of the exponent $ (<£") as a function of the eigenvalue /U. It vanishes at 
[i = 1/2 and diverges as /x — > 0. The result of the uniform approximation is plotted as a 
dashed line. 

Again we can see from the exact solution above, why the uniform approximation 
is not valid. In terms of the dominant a mode the ratio: 

Pl(0)-Pl „ J o yzH a (z)e- z ^ K 
Pk-Pl dzzH Q (z)e- z 

goes at large T towards the finite limit b a found to be bigger than one when \x < 1/2. 
One can check that setting p\ = p* k + f£ and p£(0) = p* k + bf k into the equation 
(plD, one obtains that the relaxation towards the fixed point is now like T~ x ^ with 
X a (b) = 2(6 — {b + l)/i Q ). Since b a > 1, the correct exponent <3> Q = X a (b a ) is bigger 
than the naive exponent A a = X a (b = 1) . 




5. complex \x 

Next we turn to complex fi. The matrix M being real, if it has complex eigen- 
values they will come in pairs, \i a and fi* a corresponding to complex conjugate 



17 



3.0 

2.5 - 

2.0 - 
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0.5 1.0 

FIG. 4. Plot of the exponent <1? + as a function of the eigenvalue [i. It diverges at 
/i = l/2 

eigenvectors p£ and ip%)*- In the uniform approximation they will combine as 
f£ = r-- Re ( A a) rfe cos(/m(A a ) In T + <f>%) and correspond to oscillatory (growing or 
decaying) solutions (where A = 2(1 — 2/x) is complex). This situation happens in 
cyclic reactions, examples of which will be given below. Within the uniform approx- 
imation the fixed point is stable to this perturbation if Re(fi) < 1/2 and unstable 
otherwise. This however turns out not to be correct. Indeed the correct exponent 



$ (now complex) is determined by the above equation ( 52|) associated to fj,. One 
has = <£(a0*, and to each pair of /i, //* one can associate one (or two in some 

cases) pairs of $, $* also corresponding to oscillatory solutions. The oscillation 
frequency in the lnT variable is now given by im($), and the stability being now 
determined by i?e($) (i?e($) > corresponds to a stable eigenperturbation, while 
i?e($) < correspond to an unstable one). 

Interestingly, _Re($(/i)) is a decreasing function of Im(fi) > as it is increased 
from 0. Thus the region of stability in the complex fi plane is different from the one 
inferred from the uniform approximation. It is represented in (|5|). One notes that for 
complex eigenvalues one can have a fixed point naively stable (with i?e(/x) < 1/2) 
which is in reality unstable (with i?e($ < 0). 



C. asymptotic dynamics: conclusion 

Thus we have solved the problem of the dynamics near the asymptotic states of 
the form (pTT) for arbitrary reaction diffusion process. Let us summarize the results. 



When the eigenvalues /i Q of the matrix M ( ]23( ) are real, the stability of the fixed 
point is determined by the naive argument : a fixed point is stable if fi a < 1/2 for 
all a = 2, ... The decay exponents $ a are obtained in terms of the fi a as : 

$a = F[pL a ] (54) 
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FIG. 5. Stability diagram: the solid line in the complex plane of delimitates the 
region of instability (Re(&) < to the right of the line) from the region of stability ( 
Ee($) > to the left) 



where the function F[x] is defined implicitly by the single root of the equation 
U[—2x, 2 + F[x], 2] = (here < x < 1/2) and represented in Fig. (g). 

In terms of these exponents, the large time behaviour for the concentrations of 
the species is found to be 



Pk =P k 



n r 



a=2,7V-l 



(Tlnty 



+ ... 



T 2 (lnt) 2V T(lnt) 



+ 



(55) 
(56) 

(57) 



where the 0(1/T) correction in nr comes from the convergence of the landscape 



In addition there are corrections to ( 5q) which decay must faster, exponentially 



in r (i.e algebraically in time). The b k are constants, depending on the initial 
condition. The formula (|55|) can also be used to relate two late times. If the system 
is very near the asymptotic state at t', with Pk(f) = p* k + e k , (|55|) holds at t with 



bk — e aPt(T\nt') where we recall that the p k are the eigenvectors of M and 

For practical applications it is useful to note that ratios of concentrations of 
different species involve only the exponents $ a . On the other hand, because of the 
factor nr, the relaxation of the concentration of a single specie k to its asymptotic 
form is controlled (provided pi > 0) by the exponent min(l, $) (where $ is the 
minimum of all exponents $ a appearing in the corresponding formula flSSj ) for nfc(i)). 
The formula is even more interesting in the case p* k = (i.e if the specie k disappears 
in the reaction) since then the first correction becomes the dominant decay and one 
has at large time that n k {t) ~ l/(Tlni) 2+ *. Examples of such cases are studied in 



Section VD 
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Let us stress again that the difference between the exact value $ Q and the uni- 
form approximation value A a is due to the fact that near the asymptotic states the 

p r (0) 

ratios differ from one, i.e the valleys to be decimated do not have the aver- 
age distribution of species: there is a non trivial mixing between valley heights and 
concentration of species, missed by the naive argument, and responsible for the non 
trivial relaxation exponents found here. 



When some eigenvalues fj, a of the matrix M ft2H\) are complex the fixed point is 
stable provided all /i Q lie in the part of the complex plane on the left to the curve 
represented in Fig. (|5]). The specie concentrations then relax with oscillations as: 

n k {t) = l -n T {pl+ J2 ( T J%M* a ) cos ( /m (^) MTlnt) + <t>%) + - (58) 

a=2,jV-l ^ ' 

Finally, in the case where /j, — 1/2 (or more generally on the line i?e($) = 0) 
linear analysis is insufficient to determine the evolution of the system, and one must 
study the full non linear RG equation (|15"D, which goes beyond the present study 



To conclude this Section, let us recall that the results derived above concern, 
strictly speaking, the effective dynamics described by the RG equation (|15"|). As 
was discussed in great details in Ref. |)U| for the single particule diffusion, there are 



corrections in the real dynamics, with respect to the effective dynamics. Indeed, in 
the effective dynamics the whole thermal packet jumps at hit = T over a barrier T, 
while in the real dynamics typically a fraction of a thermal packet (which can be 
written as 1 — exp(— e™ 1 ^ 1- ^)) has not yet jumped at time t. Since the distribution 
of barriers becomes broader and broader, this generates corrections which at large 
time are only subdominant for most quantities (at most 0(1/T)) coming typically 
from rare events such as degeneracy of order 0(T) between neighboring barriers. 
They become dominant only for certain quantities, such as the width of the thermal 
packet, which have vanishing leading order in the effective dynamics. In Ref. |Hj 
the corrections to first order in 0(1/T) were evaluated and found to originate from 
three rare events (a) valleys with degenerate minima, (b) almost degenerate barriers 
(c) valley just being decimated with a barrier T + e (see fig. 7 of |60||). 

A similar detailed study of the rare events in presence of reaction processes can 
be performed but goes beyond the present paper. With similar arguments as in |Kj 
we do not expect any correction to the leading order of the quantities computed in 
this paper. In principle, subdominant corrections could add to the subleading terms 
computed above, e.g. in ([55]). They are certainly at most of order 0(1 /V) (and thus 
cannot affect any decay as T - * with $ < 1) but it is likely that they are even of 
higher order. Indeed most of these corrections (e.g. (b)) come from single particle 
diffusion and can be reabsorbed into the global factor nr. Other events (such as (a)) 
cannot affect specie concentrations. Although this point deserve further study, it is 
likely that the corrections from the real dynamics to the convergence to asymptotic 
states obtained in this Section using is subdominant. 
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IV. EXAMPLES OF PROCESSES 



Up to now we have only applied the general results only to the process ([|). We 
give here several examples of other processes, starting with a process which exhibits 
a dynamical phase transition and to which the general results apply directly. Then 
we present other cases which raise interesting questions which go slightly beyond 
our present analysis. 



A. a dynamical phase transition 




Let us consider the system involving two species A and B and the empty state 

A + A-* A B + B -> B A + B ^ (59) 

The solutions of equation ([H]) for the fixed points are the empty state E = 
(p* A = 0,p* B = 0,Pq = 1), the A-phase (p* A = l,p* B = 0,pl = 0), the B-phase 

(Pa = O.Pb = i.Po = °)> and the fixed P° int C = (Pa = \,P*b = hP*o = §)■ 

The A phase and the B phase are attractive fixed points with eigenvalues ^2,3 = 

corresponding to asymptotic exponential decay in F = Tint (i.e a power law in 

time) of the other specie. The empty state is a repulsive fixed point with eigenvalues 

A*i 2,3 = 1- The critical point C is attractive for symmetric perturbation 5p A = 5p B = 

— Y corresponding to eigenvalue /i 3 = 0, but instable with eigenvalue \i2 = \ for any 

dissymetric perturbation bp a 7^ &Pb- This corresponds to the exponent $~(f) = 

—0.761258 (the other root is $ + (|) = 3.51853). Since it is globally attractive over 

the critical manifold, this fixed point controls the dynamical transition from the A 
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phase to B phase. Thus we conclude that if one starts with a system of A and B 
in almost equal concentrations, the difference ~ PB(t)\ (or equivalently the 

relative concentrations of A and B) grows with time as: 

\PA(t) ~PB(t)\ ~ (lnt)" v = 0.761258 (60) 

or, equivalently, the differences of absolute concentrations decay as \N/Sf)— N B (t)\ ~ 
(lnt)~ 2+!/ , i.e more slowly than the decay of both concentrations of A and B, which 
itself behaves as (lnt) -2 . The system eventually reaches a broken symmetry state 
where either A or B predominates after a time Ur which scales as: 

t br ~ e c(t')\PA(t')- PB (t')\-^ (61) 

where t' is a (shorter) reference time scale and c(t') a (£' dependent) constant. Note 
that the uniform approximation would predict v — 2/3 significantly smaller that 
the exact result. Finally, the asymptotic final decay of the minority phase is fast, 
exponential in T (as /x = at either A ot B fixed points). 

Restraining from giving further examples among the large number of possible 
processes with similarly interesting behaviour to which our general results readily 
apply, we now turn instead to cases where open questions remain. 



B. Reaction with a marginal fixed point 



Let us now consider the reaction: 

A + A^O B + B -> B A + B ^ (62) 

The fixed points are the empty state E{p* A = 0,p* B = 0,Pq = 1) which is unstable 
with eigenvalues /i2,3 = 1, the B phase (p* A = 0,p* B = l,p$ = 0) which is fully 
attractive with eigenvalues ^2,3 = 0, and the fixed point U (p* A — \,p* B — 0,Pq = \) 
which is attractive with /i 3 = along the axis p B = and marginal with /i 2 = 1/2 
for any perturbation 5p B > 0. Going beyond the linear approximation, we find in 
the uniform approximation that TdrP B — 2(p B ) 2 , i.e. a small initial proportion p B ° 
grows very slowly as 

and thus the time scale t eq where the proportion of B becomes finite grows like 
r = T\nt eq ~ e 1 ^ 2pB ^ for small p B °. In the real dynamics, the exponent = 1/2) 
being we also expect a kind of marginal behaviour near the fixed point U. A full 
study of this behaviour is an interesting question which goes beyond the present 
paper. 
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C. cyclic reactions and complex eigenvalues 



Let us study the reaction: 



A + A^A B + B^B C + C^C (64) 
A + B^B B + C^C C + A^A (65) 



The solutions of fl2"l~|) are the three pure phases pa — 1, Pb — 1, Pc — 1 an d the mixed 
state Pa = Pb = Pc = §■ The pure phase = 1 is stable (/i = 0) in the direction 
Spc = —Spa > 0, and unstable (/x = 1) in the direction 5ps = —8pa > 0. The 
mixed fixed point Pa = Pb = Pc = \ has complex eigenvalues /i 2 ,3 = (1 ± i/\/3)/2, 
leading to purely imaginary naive exponents A2,3 = ±i2/\/3. As can be seen on 
Fig the exact convergence exponents $2,3, solutions of fl5*2"|), have a negative real 
part, and thus the fixed point pa = Pb = Pc = | is also unstable. This shows that 
the asymptotic behaviour of the system is more complex than being described by a 
fixed point of type (pl|). 

In fact, going back to the equations (|3~2"D of the uniform approximation, and 
eliminating p c = 1 —Pa~Pb, we obtain that the flow equations for the two variables 
(Pa,Pb) take the "divergence free" form 

TdrPA = 2d PB f{p A ,PB) (66) 
Td r p A = -2d PA f(pA,PB) (67) 

with /(pa,Pb) = PaPb({—Pa—Pb)- As a consequence, all starting points where the 
three concentrations (pa,Pb,Pc) are non-zero, belong to closed flow lines p^ ni ^ orm (r) 
of constant value of f{pA,PB)- Thus in the uniform approximation the asymptotic 
behaviour is always cyclic. 

This however does not carry to the real dynamics, beyond the uniform approx- 
imation, since that one can check that these cycles p^ m ^ orm (r)e _,?1_r ' 2 are not so- 
lutions of the RG equation fll~5|) . Thus the question of determining the asymptotic 
behaviour of this problem is still open. A more complex cyclic solution, or a new 
non trivial fixed point are among the possiblities. 



We close this section by noting that one can also expect from pT3|JT^JT^| that 
reactions with a large enough number of species have chaotic solutions at the level 
of the uniform approximation. It would be interesting to investigate whether such 
chaotic solutions could also exist in the RG and in the exact dynamics of these 
disordered reaction diffusion problems. 



D. Reaction with an infinite number of states 

We now consider the much studied A + B — > inert reaction, which, in the absence 
of disorder is known to lead to segregation |3(J of the two species. 

A + A^ A + A B + B -)• B + B A + B ^ (68) 

We introduce the notations A = 0, A rn = mA and A_ m = mB (m > 1). The 
possible contents for a valley are now the A m with m e Z and thus their number is 
infinite. The reaction rules become with these notations 
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A k + A p -> A fc+P (69) 

So for the RG procedure, it is convenient to associate to each valley an auxiliary 
variable m representing the contents of the valley, and to write the RG equation 
for the probability distribution P r (zi, Z2',m) where the RG rule for the auxiliary 
variable m upon fusion of valleys simply reads m = m\ + m 2 . We can use the result 
of the Appendix of ( p4"| ) (for the same RG rule of an auxiliary variable) and obtain 
the scaling 

< m 2 >~ T 2 (70) 

Thus we find that charges of order r = Tint of both signs (i.e groups of size of 
order r = Tint of particles of the same specie) will accumulate near the bottom 
(in a packet of typical size 0(1)) of each renormalized valley. These packets will be 
separated by a large distance of order (Tint) 2 . The total number of particles will 
still decay, as l/(Tlnt). This asymptotic state thus still presents strong features 



reminiscent of the segregation observed in the pure case |30|. By contrast with the 
pure model, here several packets of A can also be found in successive neighboring 
valleys. 



V. PERSISTENCE PROPERTIES 

We now study persistence properties in the reaction diffusion models defined in 
(P~B|). As explained in the Introduction one is interested in computing probabilities 
that some type of event has not occured between time and t. The decay with time 
t of these probabilities usually involves new non trivial exponents which characterize 
the non equilibrium dynamics. Since they integrate over time the behaviour of the 
system, they are usually hard to obtain analytically, even in the pure systems. For 
reaction diffusion models in a random environment the following types of persistence 
probabilities can be defined, and will be studied in the following corresponding 
sections. 

A - The simplest persistence observable is the probability U.(t) over runs and 
environments that a given point x = has not been crossed by any particle between 
time and t in a given run. The decay of II(t): 

11(f) ~ l(t)- e (71) 

where l(t) is a characteristic length of the diffusion process, here given from (|l~0|) as 
Z(t) = |(Tlnt) 2 , defines the persistence exponent 9. This is the definition used in 
this paper, even when referring to the pure case, where l(t) ~ y/i, whereas another 
frequent definition is in terms of the power of f . Since here the diffusion is logarithmic 
we choose everywhere in this paper the more general definition ([TID both for pure 
and disordered problems. 

B - In the presence of quenched disorder one can also study the probability 
^th{t) over environments that a given point x = has not been crossed by any of 
the thermally averaged trajectories < x(t) > of the particles. Similarly the decay of 
rit/i(t) ~ l{t)~ e defines the exponent 9. One expects in general that 9 < 9 and here 
we find that these two exponents are quite different. 
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C - More generally, one can defined the probability that a given pattern present 
at time has survived up to time t. We study the example of the survival of domains 
(i.e. intervals between particles), which in the pure case was shown to lead to the 
definition of two new exponents [^6| called 5 and ip : 6 characterises the probability 
that a domain has survived up to time t without merging with other domains, and ip 
characterises the probability that a domain has survived up to time t with mergings 
with other domains. 

D - Finally we study the exponents 5 a and if) a characterising the probability that 
a particle A has survived up to time t, without any coalescence and with coalescences 
respectively. 



A. Persistence in a single run 

1. no crossing by any particles: exponent 9 

We start by computing the probability II (t) that x = has not been crossed by 
any particle up to t. We consider a rather general reaction diffusion process with 
a vacuum state (k — 0) and occupied states (with particles in them), k > 1. To 
solve this problem we can consider separately the two half spaces x > and x < 
and study the problem of a semi infinite system (x > 0) with an absorbing boundary 
at x = (absorbing for the states k > 1). 

For diffusion in a Sinai landscape in presence of an absorbing boundary at x = 
one defines a new RSRG with slightly new rules: the boundary RSRG, explained in 
detail in Ref. [[Kj. The first bond is by definition always ascending with an infinite 
barrier (and thus can never be decimated) and represents an "absorbing zone" (see 
Figure (|7|)). If the smallest barrier in the system at T is the third bond from the 
boundary or further the rules are identical to bulk RSRG. If the smallest barrier is 
the second bond, i.e the first descending bond, the procedure consists in eliminating 
the first valley (i.e the second and third bond) which is merged with the absorbing 
zone. 

Since the reaction rules of the species upon merging valleys are unaffected by 
the boundary, at a given T all renormalized valleys in the bulk are distributed 
independently with P k (771,772) which satisfies the equation (|15|). We now explicitly 
check that the first renormalized valley has also the same distribution. Indeed the 
probability Rk{r),r)') that the first renormalized valley has (77, 77', k) satisfies the RG 
equation: 

{Td v - (1 + V )d v - (1 + rf)fy - 2) R r k ( V , rj') (72) 

= P v k M) I ^^4(0 ))?2 ) + <XM)*^J(,V) (73) 
J k> 

+Wl k2 R r kM-)*r l 'P^.)-R r k (v,v') I d V2 J2 P k^V2) (74) 

J k> 

where the first term corresponds to the decimation of the second bond (which results 
in the increase of the absorbing zone) and the old second valley becomes the new 
first renormalized valley. The third and fourth terms correspond to the decimation 
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(a) 




o 



(b) 




FIG. 7. Illustration of the RG in presence of an absorbing boundary, (a) the boundary 
at site x = can be represented by setting Uo = — oo. (b) renormalized landscape, with 
the absorbing zone (see text). 

of the left bond of the second valley (the loss term must be explicitly written on the 
l.h.s. to keep the distribution R correctly normalized to unity). One can now check 
that R\(r],ri') = P£(r),r)') is a consistent solution. 

To compute the persistence exponent, we now define Vk(rj,rj') as the probability 
over all environments that the boundary at x = has never been crossed by a walker 
between time and T and the first renormalized valley has (i], rj', k) at T. It satisfies 
the RG equation: 



similar to (|72|) except for the first term which carries the restriction that the second 
bond can be decimated only if the first renormalized valley is empty (since, if it 
contains a particle, this particle gets absorbed by the wall, i.e crosses the origin). A 
consistent solution is simply 



(Td r - (1 + V )d v - (1 + - 2) V?(r], V >) 

= P v k M) J d V2 V r (0, V2 ) + Wl: iM V^.,0) * v Pl(,r)') 



(75) 



(76) 



(77) 



k' 




(78) 



with 
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Td r v r = -v r [ (j2 P k^V2)-PLo(^V2)) (79) 
Jm \ k , J 

We now use the fact that the system reaches for large V an asymptotic state 
corresponding to an attractive fixed point (pD|), and this leads to the asymptotic 
decay 

v r ~ r-t 1 -^ (80) 

Since the probability H(t) that the point x = has not been crossed by any particle 
up to time t on the infinite line is the square of the corresponding probability for 
the semi infinite problem we obtain: 

n(t) ~ v 2 v ~ i(ty e (8i) 

with l(t) = |(Tlnt) 2 and the result for the persistence exponent: 

6 = l-p* (82) 
As an example, we explicit the result for the process @ where p$ = ^-: 

e r = T ±r r (83) 

For r = where the particles A always merge and occupy all valleys, we recover 
the half-space exponent |0(r = 0) = | corresponding to the decay exponent of the 
probability of no return to the origin for a single Sinai walker obtained in p3.£)(| . 



2. number of particles absorbed by a wall: generalized persistence 

A generalization of the persistence exponent 9 can be defined for reaction- 
diffusion models on the semi-infinite line x > in the presence of an absorbing 
boundary at x — 0. There one can compute the probability Qr{n) that exactly n 
particules have been absorbed by the wall up to time t. Generalising the approach 
of (|T5D, we obtain that Qr(n) satisfies at large T: 

rdrQr(n) = (l-p*)[Q r (n - 1) - Q T (n)] n > 1 (84) 
r<9 r Q r (0) = -(l-]9*)Qr(0) (85) 

The RG equation for the generating function Qr(z) = ^2 n zn Qr(n) 

TdrQr(z) = -(1 - po)(l - z)Q(z) (86) 

thus leads to the decay Qt(z) ~ r _ ( 1_P(, ^ 1_z ^. Introducing the rescaled number of 
absorbed particules 

n 

(87) 



lnT 



and using as in [59,60] the saddle point method, we find after a Legendre transform, 
that the probability distribution Prob(g) behaves as 
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Prob(g) ~ T- 2w{9) 
with the generalized persistence exponent: 

2u(g) = (1 - pi) - g + g In (jz^) ( 89 ) 

For g = 0, one recovers the persistence exponent of the half space w(0) = 9/2, 
where 9 is given by fl82"D . u)(g) has a zero minimum at g a — (1 — po) which is thus 
the value that g takes with probability one as T — > oo 

- — — = (1 — Pq) with probability one as t — > oo (90) 



B. persistence of thermally averaged trajectories 

As was discussed in detail in Ref. ||60|| , thermally averaged trajectories of a single 
Sinai walker follow the effective dynamics which we now use to study their per- 
sistence properties. Fig. | illustrates the difference in the persistence properties 
between the single run dynamics studied in the previous Section and the effective 
dynamics of thermally averaged trajectories. Let us consider the case of a valley 
with a right bond of barrier T such that the point x = lies to the left of a valley 
bottom and is separated from it by a barrier less than T. In this will 

be crossed (several times) in a typical single run, while the thermal average (x(t)) 
will remain at the bottom of the valley until it jumps to the right over the barrier 
without crossing x = 0. Thus, as was found in |JIJ for the return to the origin of a 
single walker, the exponents 9 and 9 should generically be different. 




FIG. 8. If the point at x = (indicated by a dotted line) happens to lie in a renormal- 
ized valley as shown on the left figure, it will be crossed many times by a single particle, 
while it may not be crossed at all by the thermal average < x(t) > before it jumps over 
the barrier on the right. This is not the case in the situation shown on the right figure, 
where typically no crossing of x = occurs for a single particle. 

We now compute the probability Uth(t) that the point x = has not been 
crossed by any particle up to time t within the effective dynamics. Let us define for 
each valley the auxiliary variables m^mj equal to the total number of sites in the 
descending (mi) and ascending (m-2) bonds respectively, which have not been crossed 
by any particle between and t. We define the probabilities P/ C ( ? 7i, V2, m i, that 
a valley has a specie k, bonds 771, r] 2 and variables mi, m 2 . Consider the decimation 
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represented in Fig. ||. Let us denote the two valleys corresponding to bonds (1,2) 
and (3,4) containing respectively the species ki and k 2 - Upon decimation of bond 
2 the two valleys merge and the specie k\ jumps to the bottom of the valley (3,4) 
and thus goes over the bond (2) and (3) to react there with the specie k 2 . As a 
consequence, the auxiliary variable of the new renormalized bond F% = F\ + F% — F 
evolves with the rule: 



m' 3 = m 1 + 5 klfi m 2 + 5 ku0 m 3 (91) 

This is a particular case of the auxiliary variables studied in Appendix (|A[) with 
flfc = b k = 4,o and d k = 1. 

The final result is that the fraction of sites that have never been crossed by any 
particle in the effective dynamics decays as = ~ {h)~ 9 where the persistence expo- 
nent 9 is solution of the following equation involving the confluent hypergeometric 
functions U(a, b, z) 

9 U(-p* ,29,l) = U(- P l,29+l,l) (92) 

For the process (01) one has Pq = r/(l + r), and the resulting exponent 9(r) is 
plotted in Fig. |9|. Surprisingly we find that it is numerically extremely close for all 
r (to less than about one percent in relative difference) of one half the result [IT 



for the pure system which reads \9 pure {r) = — | + ^(arccos( v ^" ( - r ,^_ 1 - ) )) 2 an d is also 
plotted in Fig. |9[ The expansion for small r gives 

9{r) = 1 - 2r + o(r) (93) 
\e pure {r) = l--r + o{r) (94) 

and thus they are definitively different. In the case r = 1 where particles always 
annihilate, we obtain 

9(r = 1) = 0.380678.. (95) 

which may be compared with \9 pure {r = 1) = 3/8 = 0.375. The difference remains 
very small for all r, as shown in Fig. [T^. 

We have also generalized the calculation presented in this Section to compute 
the number of visits of thermally averaged trajectories of particles at a given point. 
It leads again to a multifractal spectrum of exponents. The calculation and the 
results are presented in the Appendix 0. 



C. Statistics of merging domains 



Here we define domains as intervals between particles. For the reaction (|4]), on 
which we now concentrate, when a domain dies (the two particules A meet) the two 
contiguous domains can either merge if the particles annihilate (with probability r) 
or remain separate if the particles coalesce (with probability 1 — r). To characterize 
the statictics of the coarsening of domains in the pure case (i.e Potts domains with 
q — 1 + -) Krapivsky and Ben-Nairn have introduced ||6[ the following definition. 
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FIG. 9. Plot of 9(r) as a function of r (solid line) and, for comparison, ^9 pure (r) 
(dashed line) 

They define Q m (t) as the number of domains at time t which have for ancestors 
m > 1 initial domains. This is illustrated on the figure [0]. 

The total number N(t) of domains remaining at time t (equal to the total number 
of A at time t) is simply given by the sum N(t) = X]l=iQm(i) an d decays as 
~ p* A /l(t), where l(t) ~ (Tint) 2 is the characteristic length at time t, given in fllCf). 
The fraction of initial domains which have a descendant that is still alive at t is 
given by S(t) = Y^ m m Qm(t) =< m > N(t). The decay of these quantities define 
two new independent persistence exponents S and ip- 

Qi(*)~ !(*)-* (96) 
S(t) ~ l(t)-^ (97) 

and with these exponents, it is expected that Q m (t) takes the scaling form 

Qft) = — i — Q ( (98) 

The scaling function is expected to behave for small z as Q(z) ~ z a where the 
exponent a is related to (S, ip) by the relation 5 = 2 — ip + (I — ip)a. Note that the 
inequalities Qi(t) < ^2 m Q m (t) < ^2 m mQ m {t) imply that ip < 1 < 5. Note that 
here, we have again defined the exponents with respect to the characteristic length 
l(t) at time t. Thus in the pure case our definition differs from the one of [IfJ by a 
factor of 2. 

Below, we obtain the exponents S(r) and ip(r) exactly for the process 
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FIG. 10. Plot of 59{r) = 

2@pure{ r ) ~ @(v)' the difference remains very small for all r. 
It vanishes at r = and r = 0.280701.. 

1. Exponent 5(r) for the process fy) 

For a domain to survive up to time t while keeping its variable m — 1, the 
two domains walls must avoid meeting each other up to time t, but they can meet 
other exterior domain walls, provided that upon meeting they coalesce and do not 
annihilate. 

Since the two domains walls must not meet, given the properties of the effective 
dynamics in the RG, the decay of Qi(t) is governed by the events such that at 
some r the two domains belong to two neighboring renormalized valleys. At all 
later times they will still belong to two neighboring renormalized valleys and no 
decimation of the two renormalized bonds separating the two domains can occur. 
As a consequence to compute the exponent 5, we can consider separately the two 
corresponding half spaces. 

For a given half space, we introduce the probability R r (r),rj') thaithe first bond 
has never been decimated and the valley is (n, r() and there is always one walker in 
the first valley. The RG equation for this quantity reads 



(r«9 r - (1 + v)d v - (1 + r)')d v > ~ 2)# r (r/, rf) (99) 

= R r (rj, .) V (^(0, ■) + (!- r)/$(0, .)) (100) 

+J R r (., 0) *„ (P r (., rj')) + (1 - r)Pj(., v'))) (101) 

POO 

-R T (V,V') / AfcX^CO,!*) (102) 



where the — term arises because the R r , unlike the Pj: is not normalized to 1, 
and one must count the loss associated with the left bond of the second valley. 
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FIG. 11. diffusion and merging of domains. To each domain present at time t (bottom 
of the figure) is associated the number m of ancestor domains in the initial state (defined at 
the top slice t = 0). Note that some initial domains die without producing any descendant 
domain at time t. 



Integrating over (77,77') one finds that R r = f dr)dr) f R r (r),r)') evolves with 

jpr p p poo 

= ~j dr l'R V M)-rp T A J d V R r ( V ,0)-rR r J drjP^rj) (103) 

corresponding to the three forbidden cases: decimation of the first bond, decimation 
of the second or third bond when both valleys are full and annihilation occurs. 

The exponent S will be given by the decay of the half space probability R(rj, f]') ~ 
r _<5 , since the probability associated the two sides will decay as the square of the 
probility for one side, i.e as T~ 25 ~ l r . Setting R(r),rj') = T~ 5 e~ r '~ TI ' p(r], 77'), and 
using p* + (1 - r)p* A = ^ one finds: 

= [(1 + rj)d v + (1 + r/)cV + (5 - 1 - 77 - 77')] p(7/, 77') + (104) 

I rv H 

d?7ip(?7i,0) + / dr]2p{i],i] 2 )) (105) 
^0 



1 +r 

There is a solution of the form p(?7, 77') = p(rj + 77') where ^(77) satisfies: 

= [(2 + V )d v + (S - 1 - 77)] p( V ) + ^— f drj'pir}) 

1 + r Jo 



(106) 



After derivation with respect to 77, one finds a standard hyp ergeo metric differential 
equation which admits for only solution not growing exponentially at 77 — > 00, the 
confluent hypergeometric function U(^, 2 + 5, 2 + 77). The boundary condition at 
77 = then determines the exponent S(r) as the solution of the implicit equation: 

2U'(^—, 2 + 5(r), 2) + (5(r) - 1)U(-^-, 2 + 5(r), 2) = (107) 
r + 1 r + 1 

Using the functional relation zU'(A, B,z) + (B-1- z)U(A, B, z) = -U(A - 1, B - 
1, z), the exponent S(r) is finally the solution of the equation 

[/(-_L l + 5(r),2)=0 (108) 
1 + r 
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The solution of this equation is plotted in Fig. |T|. 

In the case r = 0, we find S(r = 0) = 1 as expected. Indeed, in that case where 
particles always coalesce, domains cannot merge, and thus m = 1 is the only possible 
value : Q m (t) = y<5 m ,i and thus 5 = 1 = ip as in the pure case. For the case r = 1, 
where particles always annihilate, we find 

S(r = 1) = 2.53083.. (109) 

which is remarkably close to the numerical result obtained in ( for the Ising 
pure case : S pure (r = 1) = 2.54(4). This puzzling feature also holds for other values 
of r, as shown in Fig [12[ with less than about one percent in relative difference. 
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r 

FIG. 12. Plot of the exponent 5(r) for the process (H) in the Sinai landscape, de- 
termined by equation ( |108|) (solid line), and comparison with the numerical results of 
Krapivsky and Ben Nairn (Ref. 46) for the pure case (circles). 



2. Exponent ift(r) for the process (Qj 

To compute if) we need to obtain the scaling behaviour of the average number of 
ancestors of the domain < m >~ r 2 ^ 1- ^. However this is a priori difficult, as the 
variable m is associated to a domain which can extend over many renormalized bonds 
and is thus "non local" . However we can circumvent this difficulty by decomposing 
m upon the several renormalized bonds which make up a domain, in order to have 
a local rule under RG for an auxiliary variable b associated to bonds. Thus we write 
each variable m for a domain made out of q bonds, as the sum m = b\ + 62 + •• + b q , 
of new auxiliary variables, each associated to a bond. Since q does not grow with 
T, the scaling of < m > and < b > with F are identical. We define RG rules for 
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the local b variables as follows. We consider two neighboring valleys as in Fig ^| 
with bonds (1,2) containing specie k\ and (3,4) containing specie k 2 , respectively. 
One must think of the variable b as counting the number of ancestors associated 
to a renormalized bond and thus upon decimation of bond (2) the variables b, b' 
associated to the new bonds of barriers F\ + F% — F and become: 

b = b 1 + b 2 + b 3 and b' = 6 4 if h = (110) 
b = h and b' = b 2 + b 3 + 6 4 if ki = A and k 2 = (111) 
b = b x and b' = 6 4 if k x = A and k 2 = A (112) 

The first case where the decimated valley is empty is obvious. In the second case, 
where a particle A jumps from valley (1,2) to the empty valley (3,4), the ancestors 
of the domain to the right of A previously associated to the bonds (2), (3) and (4) 
must now all be associated to the bond (4). In the third case, where the two A 
meet, the ancestors of (2,3) disappear of the problem in all cases (i.e annihilation 
or coalescence of the A particles). 

Introducing the rescaled variables ft = b/I , where $ = 2(1 — ip), the fixed point 
RG equation for the valley distribution P k {j],r)' ', ft, ft') reads: 

= ((1 + rj)d n + (1 + 7/)fy + 2 + <f>{pd p + ft' dp, + 2))P k ( V , r/, ft, ft') (113) 

+ W Lk 2 [ Pk,(-Afti,^)* v Pk 2 (.,v',i33,fti) (H4) 
J Pi 

5{ft - (ft, + (ft 2 + ft 3 )5 klfi ))5(ft' - (ft 4 + (ft 2 + ft 3 )5 kliA 5 k2i@ )) (115) 

+ W Lk 2 [ PkM;PuP2)* v 'Pk 2 (0,;P3,Pd (H6) 
J Pi 

5{ft - (ft, + (ft 2 + ft 3 )5 ku9 5 k2 , A ))5(ft' - (ft A + (ft 2 + ft 3 )5 k2 , $ )) (117) 
We introduce the two first moments (k = and k = A) 

9k(v, rf) = — ^ ^ J dftdft'ftP k ( V , n', ft, ft') (118) 

where P*(j],ij') = e~ r?_r ''. 

Using the symmetry P k (r), r]', ft, ft') = P k (r{ , f], ft', ft), we find that they satisfy 
the closed system: 

= ((1 + r))d v + (1 + Tf)dr, -V-V'- ^)9o(v, V') (H9) 

+P*o f [9o(.,V , )+9o(;0)+9o(0,.)+9a(.,0)}+ P * f [g ( V , .) + g A ( V , .)] (120) 
Jo Jo 

= ((1 + V )d v + (1 + ri')d, nl -V-V'~ $)9a(v, V) (121) 

+P*a [ V [9A(.,0)+go(.,0) + g (0,.)]+p* [" g A (.,rf) (122) 
Jo Jo 
rv' H 

+P*a [9a(v,-) + 9o(v,-)+9o(;V)]+P*o 9a(0,.) (123) 
Jo Jo 

The exponent $ is determined by the condition that the solutions go(rj,ri') and 
gAiVi 7 ]') °f this system should be well behaved at infinity (i.e. should not be expo- 
nentially growing). 
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We found that setting 

g A (r ] , v / )=S A (z = V + 7 1 ') (124) 
9o(v, V') + 9o(v', V) = S (z = V + V') ( 12 5) 

allows to obtain the following closed system for the two functions S (z) and S A (z) 

= (2 + z)S' (z) -(z+ $)S {z) + 2p* f Z [S (.) + S A (z)\ (126) 

Jo 

0=(2 + z)S A {z)-(z + $)S A {z) + [ Z [S A (.)+p* A S (z)) (127) 

J o 

i.e, independently of the antisymmetric part of go(ViV') which we will not need. 
To decouple this system, we introduce two linear combinations S±(z) = c A S A (z) + 
CqSq(z) that satisfy closed equations 



= (2 + z)S' ± (z) -(z + <S>)S±{z) + u ± (r) f S ± {.) 

Jo 

where the eigenvalues, using p$ — r/(l + r) = 1 — p* A , are: 



(128) 



. , 1 r \/l + 6r + r 2 

u+(r) = - + ± ^—f- ^— 129 

±v ; 2 1 + r 2(r + l) v ; 

The only solutions of (|128 ) that are not exponentially growing at infinity are again 
given in terms of degenerate hypergeometric function 

S±(z) ocC/(l-i/±,3-$,2 + z) (130) 

The boundary conditions 25^.(0) = ^^-(-(O) finally give 

tf(-i/±(r),2V±(r),2) = (131) 

Since ip+( r ) < i ) -{ r )-, the growth of < m >oc T 2<yl ~'^ will be governed by ip + (r), and 
thus the final result is that the exponent if)(r) is determined by the equation 

U(— V , 2ip(r), 2) = (132) 

V 2 1 + r 2(r + l) ' ^ V ; ' ; V ; 

In particular we have the following expansion around r = 

if>(r) = 1 - -r + o(r) (133) 

For the case r = 1 where particles always annihilate, we find 

= 1) = 0.254821.. (134) 



which again is remarkably close to the numerical result obtained in for the Ising 
pure case: ip pU re( r — 1) = 0.252(2), a property which hold again for all r, and again 



to within less than about one percent in relative difference, as illustrated in Fig |13 
In the end we note that one can generalize the bound: 
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0.5 

r 



1.0 



FIG. 13. Plot of the exponent ij)(r) for the process (Q) in the Sinai landscape, de- 
termined by equation ( |132| ) (solid line), and comparison with the numerical results of 
Krapivsky and Ben Nairn (Ref. 46) for the pure case (circles). 



V'pure < %>ure 



discussed in 46 , to the disordered CclSG, cLSi 



ip(r) < 6(r) 



(135) 



(136) 



i.e that the exponent ip(r) is always bounded by the persistence exponent 9(r) of 
thermal averaged trajectories found in (|92|). This comes from the observation that 
a point that has never been crossed by any particle up to time t for the effective 
dynamics has to belong to domain that has a descendant still living at time t. Here 
the reverse inequality is clearly not true (for p* Q ^ 0) since a surviving domain may 
not contain any persistent site, as it can shift from its initial position, as shown in 



figure |14[ In particular we have found (|93|) : 



Q(r) = 1 - 2r + o(r) 



(137) 



Thus ip(r) and 8(r) differ already at first order in r. This is different from the 

2^,^l|], where it is found that 



case of the random field Ising model, studied in 

ip = 6 = (3 — v / 5)/4, as the situation depicted in figure [h| does not occur. 



D. Statistics of coalescing particles 



We now come to the study of persistence properties associated to a particle. 
Following the general framework presented in the previous section for the study of 
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FIG. 14. configuration where a surviving domain contains no persistent site, which 
accounts for the strict inequality ifi(r) < 9(r) found here for r > 0. 

domain merging statistics, we now introduce the number D n (t) of particles A at time 
t which have n particles A for ancestors in the initial condition. This is illustrated 
on the figure flToT). This will lead us to introduce two new exponents, 5 a and ipA, 
which have not been computed (nor defined) in the pure case and which we will 
compute here in the disordered model. 



The total number N(t) of particles at time t is simply given by the sum N(t) = 
^2 n D n (t) ~ (7(i)) -1 . The fraction of initial particles which have a descendant still 
alive at t is given by SU(£) = *52 n n -Dn(t) =< n > N(t). Again, the decay of Di(t), 
i.e of the number of particles which have encountered no other particles, and of 
SU(£) define two new independant exponents 




£>i(t)~7(t) 



(138) 
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s A {t)~l(t)-^ 



(139) 



Together with these exponents, it is expected that D n (t) takes the scaling form 

D ( t ) = _J V ( ^ ^ 

where the scaling function behaves for small z as T>(z) ~ z° A for small z where the 
exponent a a is related to (8a, iPa) by the relation 8a = 2 — i[>a + (1 — ^a)<7a- 

We will compute these exponents via the RSRG by two methods. The first one 
is direct, while the second one, presented at the end, will rely on results previously 
established in Section JT|. 

In the first method, we introduce an auxiliary variable n for each valley containing 
a particle A that counts the number of ancestors of this particle. We then introduce 
the probability Pa(i], rf, n) that a valley at scale V has (r/, rf) and contains a particle 
A having n ancestors in the initial condition. It satisfies the RG equation 



[ Tdr _ (i + rj) dv _ (i + v ') dv , _ 2 ] pjfa, r,', n) (141) 

= P V A ( V , ., n) V P r (0, .) + P V A (., 0, n) *„ P r (., //) (142) 

+P r (r ? , .) V P^O, - , n) + P r (., 0) *„ P r A (., r/, n) (143) 

+ (1 - r) [P^, ., .) V , n Pj- (0, ., .) + Pi(., 0, .) *„ >n Pj(., V, .)] (144) 



where P^ (rj, rf) is the probability that a valley at scale Y has (77, 77') and contains no 
particle, i.e. it satisfies fll5|). At this stage, the variable n is an integer n — 1, 2, 3, • • -, 
and the convolution on n is a discrete convolution. 

At large T, we know the fixed point Pq(t], rj') = pQe~ v ~ v and Xm^A^'V) 71 ) = 
p* Ae -v-v' where = = 1— p^. Thus setting P^ (77, r/', n) = p A e~ v ~ v ' A r (r? + ?? / , n), 
we find that the function A r (z,n) satisfies 

[Td r - (2 + z)«9, + z] A r (z, n) = 2p* f" dz'A r (z', n) + (1 - rKA r (., .) A r (., .) (145) 

Jo 

and we now will compute successively 8 A (r), V'a(^) for all r and the scaling function 
for r = 0. 

1. Exponent 8a (?) for the process (fjj 

To compute the exponent 8a giving the decay of the number of particles A that 
have only one ancestor in the initial condition, we need to solve the previous equation 
for the particular value n = 1. A r (z,n = 1) decouples from the other values of n, 
and satisfies the linear equation 

[Tdr- (2 + z)d z + z] A r (z,l) = 2p* [ dz'A T (z\l) (146) 

Jo 

Since D n=l (t) = N(t) ~ J dzze~ z A r (z,l), we set A r (^,l) = 

T 2 ~ 2Sa A(z). The function A (z) is then the solution of the differential equation 

(2 + z)A"(z) + (28 A ~ 1 - z)A'(z) - (1 - 2p*)A(z) = (147) 



(140) 
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together with the boundary condition at z = 

2A'(0) + (25 A - 2)A(0) = (148) 

We thus find A(z) oc U(l — 2po, 1 + 25 A , 2 + z), and the exponent 5 is determined 
by the implicit equation 

2U'(1 - 2p* , 1 + 26 A , 2) + (25 A - 2)17(1 - 2p* , 1 + 25 A , 2) = (149) 

Using the functional relation zU'(A, B,z) + (B — 1 — z)U(A, B, z) = -U{A -1,B- 
l,z), and pi = the exponent 5 A {r) is finally the solution of the equation 

2r 

U( — ,25 A (r),2) = (150) 

r + 1 

For the particular case r = 1 where particles always annihilate upon meeting, we 
have 8 A (r = 1) = 1 as it should since in this case the particles can have only one 
ancestor D n (t) = 5 n> iN A (t). In the limit r — > 0, where particles always coalesce 
upon meeting, we have 5a(?" = 0) — > +oo : indeed in this case at large T, all valleys 
contain a particle A (p^ = 1), and the probability to have n — 1 decays exponentially 
with T, since it requires that four consecutive bonds ( the two bonds of the valley 
and the two neighbors ) are not decimated. 



2. Exponent il) A (r) for the process ffjj 

To compute the exponent ipA, we introduce the rescaled variable v = r2 (^ A ) ■ 
The fixed point solution A(z, v) = P ^-] of the rescaled variables has to satisfy 

[(2 + z)d z -z + 2(l- il) A ){vd„ + 1)] A(z, v) (151) 

+2p* Q [ dz'A(z', v) + {l- r)p* A A(., .) * z>v A(., .) = (152) 
Jo 

In particular, using p^ = = rp* A , we find that its first moment C(z) = 

J °° dwA(z, v) satisfies the differential equation 

(2 + z)C'\z) + (2^ A - 1 - z)C'(z) + (1 - 2p* )C(z) = (153) 

with the boundary condition at z = 

2C'(0) + 2(Vu - 1)C'(0)) =0 (154) 

So finally C(z) oc Z7(— 1 + 2pQ, 1 + 2ip A , 2 + z), where the exponent i[) A is the solution 
of the implicit equation 

2U'{-1 + 2p* , 1 + 2^ A , 2) + 2{4> A - l)U{-\ + 2p* , 1 + 2^a, 2) = (155) 

Using again the functional relation zU'(A, B, z) + (B — 1 — 2) [/(A, 5, z) = —U (A — 
1,B — l,z), and Pq = the exponent ip A (r) is finally the smaller solution of the 
equation 
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tf(-_?-,2Vu(r),2) = 
1 + r 



(156) 



For the particular case r = 1 where particles always annihilate upon meeting, 
we have iPa(t = 1) = 1 = 5(r = 1) as it should, since in this case the particles can 
have only one ancestor. In the case r = 0, where particles always coalesce upon 
meeting, we have ^(r = 0) = : indeed the probability for an initial particle to 
have a descendant living at T is one, and thus S^t) is constant and not decaying. 



3. Scaling function 

The distribution T>[y) of the rescaled variable v — n can in principle 

be obtained in terms of the solution A(z, v) of (|151[) as 

V{y)= I™ d Vl H dTfr ^^T^ = [°° dzze- z A{z, v) (157) 

J0 J0 P*A Jo 

In Laplace with respect to u, we have that A r (z, q) = J °° dve~ qh 'A r (z, v) satisfies 

[(2 + z)d z - z - 2(1 - Mr))qd q ] A r (z, q) (158) 

2r f z 1 — r * 

+ / d^A{z\ q) + — - A r (, q) * z A r (., q) = (159) 

1 + r J 1 + r 

In the case r = 0, where particles always coalesce upon meeting, the number n 
of ancestors should have the same statistical properties as the length of a valley, and 
thus using the fixed point solution (??), we should have 

dve- qu AZ,U} = —\ e -V9c°mV9 (160) 



V* a sinh^ Jq 



Indeed we find: 



A r=0 {z,q) = \ e *(Wff«*»i>/5) ( 161 ) 

sinh ^fq 

is the solution of ( |158| ) for (r = 0, ipA — 0), and thus in this case the scaling function 
T>{y) reads 

P r=0 (^) = LT q X [~^\ = E + V2) 2 - l)e— 2 (W) 2 (162) 



j=-oo 



+oo 



k=—oo 



(-l) k+1 k 2 e- k - (163) 
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4- second method to compute 5a (t) and ipA{r) 



To compute Di(t), i.e the probability that a given particle A has met no other 
particles up to time t, we can consider this particle A as a tagged particle, give it a 
name, say X, and consider it as a new specie in very small concentration. It satifies 
the following reaction rule: 

A + X -> proba = 1 (164) 

and of course X + — > X and the same reactions for the A as before. We need 
only to work to linear order in px and we are back exactly in the case studied in the 
section |TTT] of the dynamics near the asymptotic state, of a new reaction diffusion 
(whose fixed point has p* x = 0). The corresponding eigenvalue of the matrix M 
introduced in (|22|) is /i = p$ = r/(l + r). Here < // < 1/2 which corresponds to an 
attractive fixed point at px = (since the X disappear) and with px ~ where 
$ is solution of ( [52]) for /j = r/(l + r). Since Di(t) ~ Px/F 2 ~ r~ 2(5A we recover the 



equation ( |150| ) 



To compute the exponent ipA (r) we need to consider similarly the reaction for 
the tagged particle X: 

A + X^X proba = l-r (165) 
A + X -> proba = r (166) 

and of course X + — > X and the same reactions for the A. In this case // = 
Po + — r ) = V(l + r ) anc ^ 1/2 < At < 1 which corresponds to an unstable 
fixed point at px = 0. One finds px ~ r~* where $ is solution $_ of (|52| ) for 
/i = 1/(1 + r). Since ^(t) ~ Px/r 2 ~ Y~ 2 ^ A we recover the equation ( |156j ) which 
determines ^(r). 



VI. DISORDER IN THE REACTION PROBABILITIES 

It is interesting to study the stability of our results to an additional quenched 
disorder in the reaction probabilities given by the matrix W (i.e spatial inhomo- 
geneities). We continue to consider only the rule that species react immediately 
when they encounter, but the analysis in fact also covers - in an effective way - the 
case where reaction rates are finite and with quenched disorder. We sketch in this 
Section a possible way of applying the present RSRG procedure to this case. 

Let us consider a model where the reaction probabilities are themselves functions 
of the position k2 (x). A simple example is to allow the parameter r to depend 
on x as r(x) in the process (f|). Let us examine what happens at a decimation time 
scale T = Tint. The particle in state k\ in the decimated valley jumps over the 
barrier to a valley containing ki- Since hi is typically at equilibrium at the bottom 
of its valley, the reaction is most likely to take place at the bottom of the valley 
within a 0(1) distance of it (since this is where all the weight of the particule ki 
is concentrated). Thus as time increases, the total number of sites in the system 
where reactions can typically occurs decays as 1/r 2 . In each renormalized valley 
at T there is typically a "finite" number (i.e not growing with T) of sites x where 
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reactions occur and thus a "finite" number of possible values W(x) (a notation for 
the set of Wt k2 )- For eacri valley these form a given set fixed in time. There are 
thus a priori two competing effects: the several values taken by W in a valley result 
in an "averaging" effect for the effective W of this valley. On the other hand the 
fact that this set is finite and fixed in time implies non trivial correlations between 
two reactions occuring at different times in the same valley. 

Here we will restrict to consider a toy model where we assign a single transition 
matrix W to each renormalized valley with probability P(W). It would be accurate 
in the case where in the initial distribution the W's are correlated over distances 
much larger than the typical thermal width of a packet ~ 0(1), but still small 
compared to T 2 . This problem can thus be treated by introducing the probability 
Pk(Vij V2, W) that a valley has rescaled barriers 771, r] 2 and an associated rate W. 
When two valley merge upon a decimation the new one simply keeps the W of the 
lowest one. One notes that the statistical independence is again preserved by RG. 
The RG equation is simply 

Td r P k r ( Vl , V2 , W) = ((1 + Vl )d Vl + (1 + rj 2 )d m ) + 2)P k r ( Vl , r? 2 , W) 
+W k k lM [Pl(vi,-,W) * m P fc r 2 (0,.) + P fc r l (.,0) * V1 Pl(.,m,W)\ (167) 

where J w Pk(ViiV2,W) = PkijIiiV^) and summation over repeated indices is im- 
plied. We also note that the distribution of W, Pv{W) = J2k I m m ^{Vu V2, W) 
is preserved by the RG rule, thus P (W) = P(W). Thus we have a "marginal" 
problem, since in this toy model P(W) does not flow by RG [J7B . 



One can now look for fixed points of this RG equation under the form: 

Pk(vi, V2, W) = e-^P k {W) (168) 
where the P k {W) must satisfy: 

P k (W) = W k klM P k2 {W) J dW'P kl (W) (169) 

In the case of the model (|j) with a distribution P(r) of r, one can show that a 
solution is: 



P k (r) = P(r)pl(r) (170) 

where Pq(t) = r/(l + r) and p\(r) = 1/(1 + r) are the equilibrium occupation 
probabilities for the problem with a uniform r, solution of the equation p* k {W) = 
W ki k2 p* k2 (W)p* ki (W). Such a simple solution holds in that case because of the form 
of the matrix M ( |2"2"| ) which is simply a projector onto the vector p* k . In general 
this does not hold and one has to solve the above equation. It is then possible in 
principle to perform, for an arbitrary P(W), the same study as the one done here, 
such as stability eigenvalues around the fixed point etc., which is left for the future. 

To summarize, the above result indicates that within the toy model and the 
effective dynamics, quenched disorder in the transition matrix W(x) will lead to a 
modification of the large time properties. These properties can be computed using 
the RG by assigning a an effective reaction probability matrix of each valley. They 
depend in a continuous way on the asymptotic distribution P(W). There is thus 
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an infinite dimensional manifold of fixed points in the RG sense, and the problem 
is marginal. Like all marginal problems, it is very sensitive to corrections which 
may make disorder marginally irrelevant or marginally relevant (or remain strictly 
marginal). The averaging effect may in the end make the disorder marginally irrel- 
evant, but to decide, within the real dynamics |76| and within a model with shorter 
disorder correlation length, whether disorder is actually marginally irrelevant, and 
how it flows, requires a more detailed study which goes beyond this work. 



VII. CONCLUSION 

In this paper we have studied the problem of various species of particles diffusing 
in presence of quenched random local bias (Sinai landscape) and reacting upon 
meeting. We have shown that the real space renormalization group method (RSRG), 
which has proved to be a powerful tool to study single particle diffusion in Sinai 
landscape ]59[ , can be extended in a simple and natural way to study a large class of 



reaction diffusion models. Since here also the physics is controlled by infinitely broad 
disorder fixed points, this method, as in the single particle problem, is expected to 
yield the exact large time behaviour. Focusing on renormalized valleys as well as 
on the particles (and species) contained in these valleys, and following the evolution 
of their distribution by decimation upon increase of the time scale, allowed us to 
obtain many new exact results on this problem. 

We have obtained a detailed description of the asymptotic states, such as the 
large time decay of the density of each specie, n k (t), and the spatial distribution 
of particles. It confirms that in d = 1 Sinai landscape the reaction is subdiffusion 
limited. The first step was to identify simple fixed points of the valley distribu- 
tion RG equation, which correspond - for a given reaction process described by a 
transition matrix - to possible asymptotic dynamical states. Each of these states 
is characterized by fixed fractions pi for each specie, the physical picture being 
the following. At time scale T = Tint the system consists of a set of successive 
renormalized valleys, which can be either empty, with probability p$ or contain a 
particle of specie k, with probability p* k . The separation between particles grows 
as the characteristic length l(t) ~ T 2 , and thus n k {t) ~ p* k / (Tint) 2 . The decay of 
concentration, when compared to the single particle diffusion length, leads to define 
universal amplitudes, which we obtained exactly. Also, from the exact statistical 
independence of the successive valley lenghts, the distribution of interval between 
particles (domains) was derived (and compared with some pure case results). 

To confirm that a given fixed point is indeed an asymptotic state, actually 
reached by the system at large time, it is necessary to study its linear stability. 
We have thus obtained analytically the spectrum of stability eigenvalues around 
any simple fixed point, as a function of the reaction transition matrix, thereby solv- 
ing the stability problem. The convergence towards these asymptotic states (i.e the 
attractive RG fixed points) was studied. The leading convergence towards these 
asymptotic states was found to be generically as \pk(t) — p* k \ ~ (Tint) - * with a 
non trivial solution of a hypergeometric equation (with in addition, an amplitude 
periodic in mint in the case of complex eigenvalues). In some cases, corresponding 
to $ = +00, the convergence is faster as a power law in time with non universal 
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exponent depending on details of the initial model. 

Eigenvalues corresponding to unstable fixed points, which were also determined, 
are of particular interest for reactions which lead to several distinct dynamical phases 
(i.e several possible asymptotic states). The transitions between different dynamical 
phases being controlled by such unstable RG fixed points, we have thus obtained 
the corresponding critical exponents. As an example, a process with a non trivial 
phase diagram was studied. 

We have also studied persistence properties associated to a given asymptotic 
state. As in pure systems, where it was originally defined, persistence can be studied 
for various types of patterns (single particles, domains etc.). Remarkably, for the 
disordered models at hand we are able to derive a much larger set of exact results 
than exists at present for the corresponding pure systems. 

We have first obtained, for a generic process, the decay exponent 9 for the proba- 
bility of no crossing of a given point by single particle trajectories. As noted in 
in a disordered system, persistence of thermal averages can be quite different from 
single particle persistence. Thus we have also computed the probability of no cross- 
ing of a given point by thermally averaged packets, which yields the decay exponent 
9. The properly generalized persistence exponents associated to n crossings have 
been defined, and computed. Next, we have characterized the statistics of domains, 
which can disappear or merge as time increases. Restricting, for simplicity, to the 
process A + A — > or A with probabilities (r, 1 — r), we have obtained exactly the 
exponents 6(r) and ip(r) characterizing the survival up to time t of a domain without 
any merging or with mergings respectively. We have also introduced new exponents 
which similarly characterize the statistics of the coalescence of particles. We have 
then computed them, namely Sa{t) and ^(r) characterizing the survival up to time 
t of a particle A without any coalescence or with coalescences respectively. 

We have found these new exponents as solutions of hypergeometric equations. 
For comparison, the only known analytical result in the pure case is for the exponent 
9 P ure(r) for the process @). A surprising outcome was that several exact exponents 
of the model with disorder were found to be numerically very close, for all values of 
r, to some exponents for the pure system, although they are associated to completely 
different diffusion length (l pure ~ Vi while 7(t) ~ (Tint) 2 ). Indeed we found that 
9{r) ps \9 pure (r) although they are definitely distinct, and furthermore that ip{r) ps 
^pureij) and 5{r) ps 5 pure {r) where ip pure ( r ), bpureij") - not known analytically - are 
extracted from the numerical simulation of |[46|| . The agreement in relative values 
is better than about one percent for all r. It may be that this observed numerical 
coincidence could be traced to the exact statistical independence of valley lengths in 
the disordered problem, while the so called "independent interval approximation", 
gives reasonable approximation in the pure case (but, surprisingly, poorer than the 
one provided by these new exponents). This however is far from an explanation and 
further investigation may be of interest. 

The effect of additional disorder in the reaction rates was also discussed. In a 
simple case it was found to be marginal, and thus yield non trivial modifications, 
continuously varying with the disorder distribution. The question of whether this 
result is stable to corrections resulting from the real dynamics or from disorder with 
shorter correlation length remains to be further investigated 

Although we have not considered explicitly branching BARW processes, with 
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Event 


Exponent 


Equation 


no crossing of by any particle 


9{r) 


6= 1/(1 + r) 


n = g In T particles absorbed at 


u{g) 


2uj = (1 + r)" 1 - g + g ln((l + r)g) 


no crossing of by thermal, aver. traj. 


B{r) 


eU(- J ^„29,l) = U(- 1 ^„26 + l,l) 


domain survival without merging 


6{r) 


tf(-l/(l + r),l + <J,2) = 


domain survival with merging 


ip(r) 


- 1^ - 4fSt . 2 ^2) = 


particle survival without coalescence 


5a{t) 


U(-2r/(l + r),2S A ,2) = 


particle survival with coalescence 


tp A (r) 


£/(-2/(l+r),2VA,2) =0 



Table I: Summary of the results obtained in the text for persistence exponents associated 
to the decay of the probability of each indicated event, in the case of the reaction diffusion 
process A + A — ► (proba r) A + A — > A (proba 1 — r). 



additional creation of particles, it is clear that for at least some of them the physics 
should not be qualitatively too different from the one obtained here. Indeed, since in 
Sinai disorder particles are essentially confined to the bottom of large renormalized 
wells, as long as the process is such that particles are not created out of the vacuum 
and that the annihilation reactions are sufficient to maintain the number of particles 
small when at local equilibrium in a well, the reaction can be treated very similarly 
via RSRG as for the model studied here. We have thus characterized a broad set of 
reaction diffusion models with disorder. 

Finally, it is worthwile to mention that we have also identified cases which clearly 
require a more complicated analysis going beyond the present paper. For instance we 
have given an example of a marginal reaction, which requires a non linear stability 
analysis. Also, we have given an example of a cyclic reaction for which all simple RG 
fixed points are shown to be unstable. The question of the determination of the true 
asymptotic states of this reaction is thus still open. Another open interesting, and 
maybe related, question is whether reactions with large enough number of species, 
which can lead to chaotic attractors in pure cases [|l|,|13|], will also lead to chaotic 
behaviour in presence of disorder. 



45 



Acknowledgements 

We thank D. S. Fisher for fruitful discussions as well as S. Fauve, U. Tauber and 
P. Chauve for helpful remarks. 



APPENDIX A: AUXILIARY VARIABLE FOR VALLEYS 

In this Appendix we study auxiliary variables (m) associated to bonds that evolve 
upon decimation as follows. Consider the decimation of bond (2) on Fig. |^: the 
two valleys corresponding to bonds (1,2) and (3,4) and containing respectively the 
species k\ and k 2 merge, and the specie k\ jumps to the bottom of valley (3, 4) and 
thus goes over the bond (2) and (3) to react there with the specie k 2 . It is thus 
natural to consider an auxiliary variable m which, upon decimation of the barrier 
F3 = F\ + F 3 — T, evolves with the general rule: 

m' 3 = d kl mx + b kl m 2 + a kl m 3 (Al) 

where the coefficients (a k , b k , d k ) depend on the specie k which diffuses upon the 
corresponding decimation. 

We now write the valley RG equation for Pj(r/,r/ ,m,m') : 

(r«9 r - (1 + r,)d v - (1 + rfjdr, - 2)P k r ( V , 77', m, m') (A2) 

= Km f (A3) 



rk 

J mi ,m'i ,m2 ,m' 2 



[ p ki(v, ;m,m'x) * v > P k2 {0, .,m 2 ,m' 2 )5(m' - a k2 m'x - b k2 m 2 - d k2 m 2 ) (A4) 
+P kl (.,0,mx,m'x) * v P k2 (.,r]',m 2 ,m')S(m - d kl nix - b kx m x - a kl m 2 )} (A5) 

Integrating this equation over the m variables one recovers of course the specie valley 



RG equation (I5|). 



We now define the first moment : 

G k (r]i,V2) = / mxP k {r]x,r]2,mx,m 2 ) (A6) 

J mi ,rri2 

GkiViiV^) = / m 2 P k (r]x,r] 2} m 1} m 2 ) (A7) 

J mi,m,2 

Since we are looking at the symmetric case, we have that G k (r)x, 772) = G\{r} 2 , r/i). 
We can thus write the following closed equation for Gl(r]x,rj 2 ) 

(Td r - (1 + Vl )d m - (1 + V2 )d m - 2)Gl( m , V2 ) (A8) 

= W k kiM [G l ki {^ .) * V2 P k2 (0,.) + a kl P kl (.,0) * m Gl 2 (., V2 ) (A9) 

+b kl Gl(0, .) * m P k2 (.,m) + 4^0,0) V P fc2 (.,r/ 2 ) (A10) 

In the asymptotic state we use the fixed point solution P^(r]x,r] 2 ) = p* k e~ Vl ~ m , 
and write G k (r)x, tj 2 ) = e~ Vl ~ m g k (rix, V2) and obtain the equation for the new function 
9k' 
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rdrg k (vi,V2) = [(1 + Vi) d m ~Vi + ( 1 + VzWm ~ V2]g k (vi, V2) (All) 

+W£ lM \piJ 9 kl (Vu ■) + P* ki a kl 9k 2 (-,V2) (A12) 
Jo Jo 
rm rm 

+Pl 2 K g kl (0, •) + P* k2 d kl g kl {., 0) (A13) 
Jo Jo 

Since the m variable is associated to bonds, it is natural to look for solutions 
where the function (771, 772) is a function of r/i alone. We thus try solutions of the 
form gl(r]i, r/ 2 ) = g k (r)x), where the exponent i[) characterizes the scaling of the m 
variable m ~ . 

For this to work we obtain, in terms of the matrix M defined in (0), the necessary 
condition: 

9k(v) = Wl ka pl 2 g kl {r,) = M kM g kl (v) (A14) 
together with the differential equation for g k (f]): 

rv 

= [(1+ v)d v -V- ^9k(v) + (M k , kl d kl + W^ lM p* k2 a k2 ) / g kl {.) + M kM b kl g kl (0)v (A15) 

Jo 

One can then try ^(77) = p k ip(r]), which automatically satisfies the necessary 
condition ( |A14j ) above (since pi is by construction an eigenvector of the M matrix 
of eigenvalue 1), and then the second equation gives the conditions involving two 
numbers Ai i2 

M kM pt 1 {d kl + a kl )=\ 1 p* k (A16) 
M kM p* kl b kl = X 2 p* k (A17) 

together with the differential equation for g(r]): 

= [(1 + V )d v - 77 - <%(77) + Ai il>(.) + A 2 ^(0)r/ (A18) 

Jo 

We now give two applications of this general analysis. 



1. Persistence exponent 9 

We now study the case a k = b k = 5 k $ and d k = 1 corresponding to the auxiliary 
variable ( |A1| ) needed to compute the persistence exponent 9 The conditions ( |A16| ) 
become: 

p* k + M k , p* = X lP * k (A19) 
M k , oP *o = \2P* k (A20) 



Since the rates involving the empty state (0) are given by definition by Wf = 5, 
we have M kt0 = W^p* = S kt ip* = p* k . The conditions above are thus satisfied with 
Ai = 1 +p* and A 2 = Pq and thus the problem reduces to studying the the equation 
([AT8D for g( V ): 
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o = [(i + v )d v -n- if)}g( v ) + (i + Pi) I g(.) + P* g(0)vi (A2i) 







Differentiating with respect to rj one gets: 



= (1 + v)d 2 v g(v) + (1 - - 77)^(77) + p* (g( V ) + (?(0)) (A22) 

with the boundary condition g'(0) = ipg(0). The solution of this confluent hy- 
pergeometric equation that is well behaved solution at infinity (i.e not growing 
exponentially) reads: 

g( V ) = p(0)(2 y( " P5 ' 2 "^' 1 V ) - 1) (A23) 

The boundary condition at rj = then leads to the following equation for the 
exponent if) governing the scaling of the m variable m ~ T^, as a function of p* : 

U'(-p* , 2 - V, 1) = |c/(-Po, 2 - V, 1) (A24) 

Using functional relations of the confluent hypergeometric functions {7, we finally 
obtain that the fraction of sites that have never been crossed by any particle in the 
effective dynamics decays as = ~ r^~ 2 ~ (lr) d where the persistence exponent 

9 = (2 — if))/ 2 is solution of the equation 

U(-p* , 29, 1) = U(-p*, 29 + 1, 1) (A25) 



2. Number of thermal packets seen by a given point 

We introduce the bond variable m(n) which is the number of sites on the bond 
which have been crossed exactly n times by a particle (any non empty state) in 
the effective dynamics (i.e by a thermally averaged trajectory). It satifies upon 
decimation of bond (2) with the same conventions as above: 

m' 3 (n) = rrii(n) + m^in — 1) + 777,3(72 — 1) if k\ 7^ (A26) 
777/3(77.) = 7771(77) + 7772(77) + 7773(77) if k± = (A27) 

Introduces the generating function m(z) = Y2n=o°° ^{n)^-, the rule becomes: 

777/3(2:) = d kl (z)mi(z) + b kl (z)m 2 {z) + a kl {z)m 3 (z) (A28) 

which, for fixed z, is the same rule as above with now a k (z) = b k {z) = S k) o+z{l—8 kj o) 
and d k = 1. The conditions ( |A16| ) become 



(1 + z)p* k + (1 - z)M kfi p* Q = X x pl (A29) 
zp* k + (1 - z)M k , p* = X 2 p* k (A30) 

and thus using again W^ = 6 k) i, these conditions are satisfied with \i(z) = 1 + z + 
(1 — z)pl and A 2 = z + (1 — z)p^., i.e. we only need to perform the substitution 
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pi — > z +Po(l — z) in the previous solution to obtain the equation for the exponent 
9(z) governing the scaling of the ratio ~ (lr) ■ 

9(z) U(~z - p (l - z), 29(z), 1) = U(-z - p (l ~z),l + 29(z), 1) (A31) 

The probability that a given point has been visited by n thermally averaged 
trajectories up to time t, is thus obtained in the rescaled variable g = -^f as: 

Prob(g) ~ (Zr")^ (9) (A32) 
It decays with the exponent 9(g) obtained through the Legendre transform 

29(g) =29(z*(g))+gln(z*(g)) (A33) 

where z*(g) is the solution of + - = 0. 

One can compute simply the value g a that g takes with probability one as T — > 
oo. It is given by 9(g a ) = = ^\ g=ga . This gives g a = -2^\ z=1 , and thus 
differentiating ( |A31|) with respect to z and taking z = 1 we finally get 

<" - (1 " ■O gKW-'i-ui) = 1^ " p;) (A34) 

where we have used the notations Ui(a,b,z) = d a U(a,b,z) and U 2 (a,b,z) = 
d b U(a,b,z). 



APPENDIX B: THE PARTICULAR CASE OF "ASSOCIATIVE 

PROCESSES" 

It turns out to be useful to introduce the notion of "associative processes" : these 
are processes such that the outcome of a sequence of reactions does not depend on 
the order in which it was performed, i.e. such that the rates W satisfy: 

pia rr p 2 p3 rr P2P PlP3 K^^J 

for all k,pi,p2,P3 (contraction over a and (3 is implied). This means that the prob- 
ability of (p2-Pz)-P\ — k is identical to the probability of (pi-Pz)-P2 — k (a.b denotes 
the result of the reaction of a and b). For example, the process defined in (f|) is 
associative. 

An important property of associative processes is that their matrix Mp2| satisfies: 

M 2 = M (B2) 

and thus the eigenvalues \i a have only two possible values : or 1. 

For the RG, these processes have also the following interesting property : the 
subspace of valley distributions of the form 

P£(Vu m) = W r k ur2 Pr( Vl )pr( m ) (B3) 

is conserved by the RG (|15|), provided that the bond distribution P[(v) satisfies the 
bond RG equation 
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TdrP^v) = ((1 + V)d v + l)P£(v) + W^ k3 W^ k2 Pl(0)Pl(.) * v P£(.) (B4) 

The bond RG equation (|B4|) can in fact be interpreted to characterize the fol- 
lowing modified reaction diffusion process, that we call "the bond-reaction diffusion 
process" : one associates to the bottom of each bond (i.e the point of lowest energy) 
a specie in one of the possible "states" and define probability distributions Pk{z) for 
the bonds. We consider two consecutive valleys made with the bonds (1,2) and (3,4). 
Initially the bonds (1,2,3,4) respectively contain the species ki, k 2 , k 3 , k 4 . Upon dec- 
imation of bond (2), the bond diffusion process is defined as follows in three steps 
(i) First the two species hi, ki on bonds (1,2) react to give another state k' with the 
rates W ki k2 (ii) the new specie k' diffuses towards the bottom of the bond 3. (iii) 
the species (A/, k%) react at the bottom of bond 3 to give a new species k' 3 with the 

k' 

rates W k , 3 k ). For comparison, it is useful to recall the corresponding real dynamics 
with valleys : initially, the valley (1,2) contains some species k', the valley (3,4) 
some species k". Upon decimation of bond (2), the species k' diffuses towards the 
bottom of the valley (3, 4) and reacts there with k" to give k with probability W k , k „. 
Thus, in the end, the physical content (the specie) of the renormalized valley for the 
bond-diffusion process is k with probability W k , k W k ? k3 i whereas in the original val- 
ley process, the final result is k with probability W$ k „W k3 k4 . The two descriptions 
are thus equivalent in that sense only if the rates satisfy the associativity condition 
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